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O ! Abstract 

Considering quantum cosmological minisuperspace models with positive 
potential, we present evidence that 

(i) despite common belief there are perspectives for defining a unique, 
naturally preferred decomposition of the space JH of wave functions into two 
subspaces iff ± that generalizes the concept of positive and negative frequency, 
and that 

(ii) an underlying unitary evolution within these two subspaces exists and 
may be described in analogy to the representation of a geometric object in 
local coordinates: it is associated with the choice of a congruence of classical 
trajectories endowed with a suitable weight (such a setting is called WKB- 
branch). 

The transformation properties of various quantities under a variation of 
the WKB-branch provide the tool for defining the decomposition. The con- 
struction leads to formal series whose actual convergence seems to require 
additional conditions on the model (related to global geometric issues and 
possibly to analyticity) . It is speculated that this approach might relate to 
the refined algebraic quantization program. 

PACS-number(s): 98.80.Hw, 04.62. +v . 
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1 Introduction 



Minisuperspace quantum cosmology bears some formal similarity to the quantiza- 
tion of a scalar particle on a curved space-time background. The Wheeler-DeWitt 
equation in the former corresponds to the Klein-Gordon equation in the latter type of 
models. Therefore, concepts, methods and experiences that have originally emerged 
in particles quantization issues are often taken over to quantum cosmology. One of 
the most prominent statements that have a clear meaning in particle quantization is 
that in the generic case, i.e. if the wave equation does not admit any symmetries, the 
is no unique decomposition of wave functions into positive and negative frequency 
subspaces. This statement has been transferred to quantum cosmology although in 
this field it is much less clear what positive and negative frequencies are and what 
other concepts might serve as generalizations thereof. Also, it is often believed that 
the space of wave functions cannot — in a natural way — be decomposed into a 
direct sum such that the Klein-Gordon scalar product is positive definite on one 
and negative definite on the other subspace. However, strictly speaking, this type 
of belief is true only if one refers to the construction of a decomposition by means 
of underlying local symmetries. 

The purpose of this article is to show that (assuming the potential in the wave 
equation to be positive) there are still perspectives in a number of models for the 
construction of a preferred decomposition generalizing the concept of positive and 
negative frequency. In any model admitting the application of the procedure pro- 
posed here, this decomposition is not based on a local symmetry but rather on global 
issues. It gives rise to a Hilbert space structure at the fundamental level of solu- 
tions to the (minisuperspace) Wheeler-DeWitt equation. Furthermore, in situations 
where such a preferred decomposition exists, a unitary evolution within the two sub- 
spaces is exhibited in a local context, depending on the choice of a classical action S 
and a weight function D. Calling the pair (S,D) a " WKB-branch" , the hermitean 
evolution operator emerges in some analogy to the way tensor components appear 
in coordinate systems, namely as representatives of some mathematical "object". 
The proposal we formulate is understood as a general strategy whose realization 
in concrete models will face specific problems that might be easy or hard (or even 
impossible) to solve. 

A further goal of this article is to speculate on the relation between the pre- 
ferred decomposition and structures that appear in the context of other quantization 
methods, in particular the refined algebraic quantization. This method incorporates 
global issues from the outset and ends up with a positive definite inner product on 
the space of physical states. Uncovering a relation to the more conventional quanti- 
zation in terms of the Klein-Gordon scalar product could eventually lead to a more 
fundamental and unified understanding of what happens when one "quantizes" a 
theory. 

The paper is organized as follows: In Section 2, we write down the wave equation 
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and the two scalar products (Klein-Gordon and L 2 ) at hand. Also, we briefly review 
how congruences of classical trajectories in minisuperspace are described by the 
Hamilton- Jacobi equation and its solutions. After a discussion of the decomposition 
issue and the refined algebraic quantization for the case of the Klein Gordon equation 
in flat space-time, we formulate the general problem. 

Section 3 is devoted to the detailed presentation of a formalism serving to analyze 
the wave equation in a novel way. Any WKB-branch may be used as a "background 
structure" with respect to which many decompositions are possible, each one as- 
sociated with a solution of a differential equation for an operator called H. This 
differential equation basically states that in each of the two subspaces the wave 
equation is of the Schrodinger equation type, H serving as evolution operator. The 
key idea of this procedure is to generalize the notion of outgoing and incoming orien- 
tation of classical trajectories to the quantum case. At first sight, this construction 
is not unique, and we end up with a variety of ways to do so. 

In Section 4, we analyze the transformation properties of various quantities when 
the WKB-branch used as "background structure" is changed infinitesimally. Sec- 
tion 5 is devoted to the formulation of our proposal. At a formal level we show 
that — within each WKB-branch — one may choose a solution H of the above- 
mentioned differential equation such that all decompositions emerging in different 
WKB-branches agree with one another. This opens the perspective to unambigously 
single out one preferred decomposition of the space of wave functions. The question 
whether the formal (power series) solution for H actually convergences (or may con- 
sistently be regularized) is addressed to the concrete models upon which the strategy 
shall be applied. We explicitly display the first few terms in the expansion of H (and 
some redefined version /C which makes explicit the unitarity of the evolution). In 
semiclassical contexts, these terms may be expected to approximate the actual solu- 
tion and should thus unambigously describe the decomposition and the underlying 
unitarity. 

In Section 6, we comment on the structure exhibited by the formal solution of 
the decomposition problem and give some speculations on its nature. It is clear that 
- whenever it exists — the solution refers to the global structure of minisuperspace 
(possibly by means of some average over all possible decompositions), and there are 
hints that some analyticity structure might play a role as well. In Section 7, we dis- 
play simple examples for which the operator H may be expressed in terms of known 
functions. Furthermore, using a toy-model situation, we show that the obvious 
heuristic argument against a preferred decomposition (the difference of asymptotic 
incoming and outgoing decompositions on a time-dependent background) does not 
necessarily apply. Section 8 is devoted to speculations about relations of our por- 
posal to the refined algebraic quantization scheme. We give a heuristic argument 
(whose basic ingredients date back to the literature of the Sixties) that the mere co- 
existence of the positive definite physical inner product constructed in this approach 
with the indefinite Klein-Gordon scalar product on the space of wave functions sin- 
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gles out a preferred decomposition. It is however not clear whether it coincides with 
the decomposition emerging from our proposal or whether it is yet another "natu- 
ral" solution of the decomposition problem. Finally, in Section 9, we comment on 
possible generalizations to the case of not necessarily positive potential and point 
out that this issue might touch upon the significance of analytic continuation and 
the role of the so-called Euclidean domains in minisuperspace. 

Summarizing, our approach intends to clarify the relation between structures 
emerging in different quantization schemes. At a formal level, we can offer a proce- 
dure that starts from the conventional setting of the Klein-Gordon scalar product on 
the space solutions to the wave equation but transcends this framework and exhibits 
a Hilbert space structure at the fundamental level, without making use of semiclas- 
sical techniques and other approximations. Although constructed in a local context, 
this structure is likely to be based on global issues, and it might be related (or even 
identical) to the physical Hilbert space as constructed in the more advanced refined 
algebraic quantization method. 

In order to improve the readability of this article, we list the most important 
notations used: 



(S, D) Classical action and wave function prefactor, defined in (|2.5| ) 

and Q3.6| ). The pair (S,D) is called WKB-branch. 

t —S, classical evolution parameter along trajectories, serving 

coordinate in minisuperspace Ai. 

Q Indefinite Klein-Gordon type scalar product defined in (|2.8| ). 

( | ) Positive definite scalar product, defined by ( P-13|) in terms of 

a WKB-branch. 

M Suitably defined vector space of solutions to the wave equation 

fl2.2|), such that Q exists on M. 
M ± Subspaces of M, generalizing the notion of positive /negative 

frequency. 

h Operator characterizing the wave equation in the form ( |3.9| ). 

defined in CTHCT- 



H Operator in the evolution equation Q3.23D for wave functions 



in M + , subject to the system of equations (|3.18|) - (|3.20|) . 
/C Hermitean evolution operator in the Schrodinger type equation 

( |3.38|) for wave functions in M + , defined in (|3.40| ). 



e Formal book-keeping parameter which can be inserted into most 

equations of this paper; its actual value is 1, and it is introduced 
in Section 5. 

7~C phya Physical Hilbert space as constructed by the refined algebraic 

quantization method. 

Positive definite inner product as constructed by the refined 
algebraic quantization method. 
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2 Klein-Gordon type wave equation with positive 
potential — the problem 



The basic mathematical ingredients for quantum cosmological minisuperspace mod- 
els QXJ as well as for scalar particle quantization in a background space-time are 
an n-dimensional different iable manifold M. together with a metric ds 2 of pseudo- 
Riemannian signature (—,+,+,...,-1-) and a real- valued function U on M.. In many 
practical cases, these objects will be C°° or even (real) analytic, and we postpone 
further specification of their differentiability properties to speculations in Sections 
5 and 6. In local coordinates y = (y a ) we can write 

ds 2 = g a p{y)dy a dy f3 . (2.1) 

In quantum cosmological models, M. (the minisuperspace) consists of the spatial 
degrees of freedom of geometry and matter, ds 2 is the DeWitt metric and U the 
potential, stemming from the spatial curvature as well as the matter couplings. In 
the scenarios of scalar particle quantization in a background, M. is the space-time 
manifold and ds 2 its metric. The massive (Klein-Gordon) case is given by constant 
potential U = m 2 , but more general potentials are admissible as well. 

We require in addition that there exists a foliation of M. by spacelike hypersur- 
faces. (The words "spacelike" and "timelike" shall refer to the "causal structure" 
provided by the metric ds 2 . In quantum cosmology, this is of course not the causal 
structure of space-time). As will be pointed out below in this Section, it will be 
necessary in some models to impose asymptotic conditions, e.g. to choose a folia- 
tion consisting of Cauchy hypersurfaces. In any case, we can assign each local light 
cone a future and a past orientation in a globally consistent way. Any smooth ori- 
ented timelike curve in A4 is either future- or past-directed. In the case of quantum 
cosmology, one could call these two orientations "outgoing" and "incoming", re- 
spectively. The hypersurfaces providing the foliation can be labelled as T\ (A G M), 
with future (outgoing) orientation being equivalent to (or, likewise, being defined 
by) increasing A, past (incoming) orientation by decreasing A. 

As an important restriction for what follows we require that the potential U 
is everywhere positive on M. In quantum cosmology this implies that any clas- 
sical trajectory (representing a complete space-time with matter configuration) is 
everywhere timelike, i.e. either always outgoing or always incoming. Returning tra- 
jectories and turning points are excluded, as we shall see shortly. For a particle in a 
background this condition guarantees that the classical motion (i.e. a world-line in 
the space-time manifold A4) respects the causal structure provided by the metric in 
that the four-velocity is always timelike (for generality, we admit both future- and 
past-oriented trajectories). In the massive case U = m 2 , this condition is obviously 
satisfied for real m. In Section 9 we will speculate about the perspectives to relax 
the condition U > 0. 
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Denoting by V a the covariant derivative with respect to the metric, the basic 
wave equation is given by 

(- V a V a + u)^j = 0, (2.2) 

its solutions ip = ip(y) being denoted as wave functions, or "wave functions of the 
universe" . They may be considered as representations of the quantum states of the 
system. In quantum cosmology ( |2.2| ) is the (minisuperspace version of the) Wheeler- 
DeWitt equation, whereas in the particle quantization models it may be called the 
Klein-Gordon equation in a background space-time with external potential. 

We should now make more explicit what we consider to be the underlying clas- 
sical dynamics. A classical evolution of the system under consideration is provided 
by an oriented trajectory y(t) in M. satisfying the equations of motion generated by 
the action 

S=\j dt(^g aP rf-NU). (2.3) 

The "lapse function" N serves as a Lagrangian multiplier and thus does not have a 
dynamical evolution equation. Its appearance in the action expresses reparametriza- 
tion invariance, i.e. it reflects the fact that the evolution parameter t has no physical 
significance (except for its orientation, which is achieved by restricting N to be pos- 
itive). After variation with respect to the variables (N,y a ), it may be chosen as a 
function iV = N(y), thereby fixing the evolution parameter t. We will not write 
down the equations of motion explicitly but just note that the momenta p a are 
related to the "velocities" by 



l dy 



N dt 



9 a % , (2.4) 



and variation of S with respect to N yields the constraint p a p a — — U, whose quan- 
tum version is the wave equation ( |2.2|) . Sometimes, for n > 3, when passing from 
the classical to the quantum theory, the original classical potential U cl is modified 
by adding a curvature term, U = U cl + £R[g], where £ = j(n — 2)/(n — 1) and R[g] 
is the Ricci-scalar of the DeWitt metric |]]. This makes the wave equation invari- 
ant under redefinitions of the lapse iV = Q 2 (y)N, which amounts to a conformal 
transformation ~g a p = Vt 2 g a p, a rescaling of the classical potential JJ C = VL~ 2 U cl and 
a rescaling of the wave function ip = Q( n ~ 2 ^ 2 ip. In such a case the quantity U in 
( |2.3|) shall denote the modified potential, as it appears in the wave equation, and 
when talking about classical trajectories, we understand the dynamics generated by 
the action fl2.3|) with precisely the same potential U. In other words, classical and 
quantum theory are both formulated in terms of the same fundamental objects ds 2 
and U, the latter everywhere positive. In some cases with not everywhere positive 
potential, it is possible to achieve U = U cl + £R[g] > after a transformation of 
the type given above. When trying to find an appropriate conformal factor Q 2 , one 
has to make use of the conformal transformation properties of the Ricci-scalar. One 
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should however agree on a definite choice of U to work with, because the formalism 
to be developed in this paper is not conformally invariant in an obvious way. (This 
would be an interesting issue to pursue, but we will not raise it here). 

In what follows we will not make use of single classical trajectories but of con- 
gruences generated by solutions of the Hamilton- Jacobi equation. The latter reads 

Sa s a = -U, (2.5) 

whereby we use the abbreviation S a = V a S, S a = V a S. Once having chosen a real 
solution S(y) (which in general may exist only in a region Q of A4), one may identify 
its gradient S a with the momentum p a of a family of oriented curves. Having fixed 
N to be some positive function on Ai, (|2.4|) plays the role of the differential equation 



4 ^ = S a (2.6) 
N dt y ' 

for y a (t), the general solution being a congruence (i.e. a non-intersecting (n — 1)- 
parameter family) of solutions to the classical equations of motion. Redefining 
N just rescales the evolution parameter t along the trajectories (in an orientation- 
preserving way). Choosing a particular function N eliminates the reparametrization 
invariance, and thus has the status of fixing a gauge. A Hamilton- Jacobi function S 
is sometimes called " classical action" . This is because S agrees with ( [2.3|) when the 
integration is carried out along the according trajectories (and all lower integration 
bounds lie on a fixed hypersurface of constant S). 

The tangent vectors of all trajectories in the congruence are, by construction, 
orthogonal to the hypersurfaces S = const. Due to the positivity of the potential, 
the trajectories are timelike, the hypersurfaces of constant S are spacelike. The 
trajectories belonging to some S are either all outgoing or all incoming. Thus, any 
S may be associated with one of these two possible classical modes of orientation. 
Accordingly, S' = —S generates the identical congruence, but with all trajectories 
reversed in orientation. Whenever different trajectories intersect each other (e.g. on 
a caustic), or whenever the congruence contains a sequence approaching a lightlike 
curve, S develops a singularity. (An example of this latter type is provided by 
ds 2 = —dr 2 +dx 2 , U = 1 and S = \/t 2 — x 2 ). This causes Hamilton- Jacobi functions 
S in general to be defined only in some region Q of M. (although the analytic 
continuation of S across the boundary of Q plays a role in certain approaches to 
quantum cosmology). 

After having summarized how the classical evolution emerges from the structure 
provided by ds 2 and U, we return to the issue of the quantum problem. There are 
two natural candidates for a scalar product, namely 

q{^2)= [ (2.7) 
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and 

Wi, ik) = - \ [ dY, a (i/il Va V2) , (2-8) 
2 js 

where dfi = d n y \/— 5 is the covariant volume element on M. , E a spacelike hyper- 
surface and cffi a the covariant hypersurface element on E. Both scalar products 
have been considered in attempts to quantize the system and to interpret solutions 
of the Wheeler-DeWitt equation in terms of observational physics. The first one, 
q, is sometimes modified into an integration over a bounded domain of M. and in 
this form intended to be the basis for computing relative (conditional) probabilities 
for "finding" the configuration variables in the respective domain. It is associated 
with what is sometimes called "naive interpretation" of quantum cosmology and 
has been advocated by Hawking and Page 0. It is also the starting point for the 
"refined algebraic quantization" program that is used in the more recent versions 
of the Ashtekar approach [[J and about which we will say a few words below. The 
second (Klein-Gordon type) indefinite expression Q is especially designed for the 
case that ip\ and ip2 satisfy the wave equation ( |2.2| ) (in which case its integrand is a 
conserved current) and is the starting point for another range of approaches to quan- 
tum cosmology, such as the Bohm interpretation [Q and the semiclassical (WKB) 
interpretation as proposed, among others, by Vilenkin ||. By replacing ipl — > ipi 
and omitting the prefactor i in ( ^.8| ), it is sometimes viewed as a symplectic struc- 
ture on the space of real wave functions. Some results have been achieved on the 
relation between the two scalar products q and Q [§, || [/J, but they either involve 
approximations or refer to very specifiy situation and are thus not very conclusive 
at the fundamental level. 

The formalism to be developed in this article starts from the "conservative" 
point of view that some appropriately defined space M of solutions of the wave 
equation ( |2.2| ) is endowed with the structure Q, and that the physical significance 
of a wave function should be expressed in terms thereof. One could call this point 
of view "Klein-Gordon quantization" (cf. Ref. ||). When proceeding, non-local 
and non-causal issues will appear, and we will frequently speculate on relations of 
our formalism to schemes that are based on q, in particular the refined algebraic 
quantization that was mentioned already above. This method has been developed 
during the past few years by several authors (see Refs. |E|, [7|, |9|, ITOfl ; also see Refs. 



n], [L2|, [L3|, [L4| for predecessors). In our simple case, the idea is essentially to utilize 
the "delta function" 5{C) = (2n)^ 1 j d\ exp(iAC) where C is the wave operator in 
( |2.2| ), the integration actually being related to a group average. Physical states are 
defined as i/i phys = 5(C)ip, where ip is to some extent arbitrary (in particular, it should 
not satisfy the wave equation). The inner product is given by {ijj hys \(f> ph B ) ph s = 
J dfiip*5(C)<f). This naive idea can be made precise in terms of distributions within 
the framework provided by the "kinematic" (auxiliary) Hilbert space L 2 (Ai,dfi) 
(see also Ref. fll|). The main condition for the procedure to work is that the wave 
operator C is self-adjoint with respect to this structure. 

In contrast to q, the Klein-Gordon type scalar product Q is adequate not only for 
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quantum cosmology but for the scenario of a scalar field in a space-time background 
as well, since it does not need information about ds 2 and U residing in the future of 
S. It is thus compatible with the idea of a field ip propagating causally in M.. On 
the other hand, such a feature is not at all necessary in quantum cosmology, and we 
will transcend it later. For the moment, however, our procedure may be interpreted 
in terms of quantum cosmology as well as in terms of a scalar field in a background. 

There is a subtlety related to the question which hypersurfaces £ may be used 
in ( |2.8|) . When both ipi and ^2 satisfy the wave equation ( |2.2|) , the integrand is a 
conserved current. Hence, whenever at least one of the two wave functions vanishes 
sufncently fast in the asymptotic boundary region of two spacelike hypersurfaces S 
and the complex number Q(ipi,ip 2 ) is independent of whether it is computed 
using £ or However, this is a rather vague statement, since we have not pre- 
supposed a particular asymptotic structure. One may, for example, think of E to 
be a Cauchy hypersurface, and £' to be a spacelike non-Cauchy hypersurface (if 
ds 2 = —dr 2 + dx 2 , it could be of the type r 2 — x 2 = const). Whether two hypersur- 
faces differing so much in their global structure give rise to a unique scalar product 
Q will highly depend on the model — i.e. on ds 2 , U) — as well as on the space 
of solutions to fl2.2| ) admitted. The minimal requirement is that the scalar product 
Q is compatible with £ being any of the hypersurfaces T\ defining the foliation 
of A4. In some models it will thus be necessary to impose additional asymptotic 
conditions, e.g. to require all T\ to be Cauchy hypersurfaces. We assume that a 
class of " admissible hypersurfaces" has been specified so as to fix Q unambigously. 
Whenever we use a spacelike hypersurface £ to compute Q, it shall belong to this 
class. 

In the case of the Klein-Gordon equation in a flat background (ds 2 = —d^ + dx 2 , 
U = m 2 ), the indefinite nature of Q gives rise to a decomposition of the space of wave 
functions M into positive and negative frequency subspaces (particle and antiparti- 
cle modes), M = 1H + ®1H~. The spaces M + and M~ are orthogonal to each other, 
and Q is positive (negative) definite on M + (1H~). Complex conjugation maps 
into each other. In principle, there are many such decompositions, related to each 
other by Bogoljubov transformations. However, in flat space-time the existence of 
a timelike Killing vector field as well as the action of the Lorentz group ensures 
the standard (Lorentz invariant) decomposition to be singled out as preferred one. 
Moreover, when a Lorentzian coordinate system has been fixed, the wave equation 
reduces to a Schrodinger type equation with unitary evolution within each subspace. 
One thus encounters two Hilbert space structures (M + ,Q) and (JH~,—Q) within 
each of which an evolution of the "ordinary" quantum mechanics type (although 
with non-local Hamiltonian) takes place. Fixing some Lorentzian coordinate system 
and writing the wave function as ip(r,x) = x{ T i x)e~ tmT , the evolution in 1H + (or 
likewise the definition of M + ) is given by id T x = Ex with E = [m 2 — A) 1 / 2 — m. 
This operator may be well-defined in the momentum representation (i.e. by means 
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of the Fourier transform). The Schrodinger type wave function \ appears as pref- 
actor of e tS with classical action S = —itlt, which in turn is associated with a 
congruence of classical particles at rest in the Lorentzian coordinate system chosen. 
Since this coordinate system was arbitrary, any congruence of particles in parallel 
uniform motion provides a possible "background" with respect to which a unitary 
evolution equation may be written down. It should be remarked that this reasoning 
corresponds to a "one-particle interpretation" and is not the modern way to pro- 
ceed when treating the quantized scalar particle. However, in quantum cosmology 
the wave function if> is usually not promoted into an operator — as long as third 
quantization is not envisaged — and one may look at it from the point of view of a 
" one- universe interpretation". We will, in this paper, try to generalize the pattern 
provided by unitarity and decomposition to non-trivial backgrounds (A4, ds 2 , U) as 
far as possible. If one likes in addition the space M of wave functions to form a 
Hilbert space, one may reverse the sign of Q in M~ and obtain an overall positive 
definite scalar product. Although such a procedure might seem somewhat arbitrary, 
it is completely compatible with the Lorentz symmetry. Moreover, the refined alge- 
braic quantization approach [K]], when applied to the Klein-Gordon equation in 



Minkowski space @ as well as in a flat space-time with toroidal spatial sections || 
(i.e. M. = 1R x S 1 x S 1 x S 1 with arbitrary radii for the S 1 -factors) precisely yields 
the structure (M + , Q) © (1H~, —Q) as the Hilbert space (7i phys , ( | ) phys ) of physical 
states. 

A typical way how the conceptual problems in quantum cosmology arise is that 
a generic background (ds 2 , U) does not admit a preferred frequency decomposition 
of wave functions in an obvious way. This problem persists even when the potential 
U is positive, and the notion of outgoing and incoming trajectories (i.e. the classical 
analogue of mode decomposition) is unambigously well-defined. Also, the split of 
the wave equation into two (mutually complex conjugate) unitary evolution equa- 
tions is usually achieved only in a semiclassical context H where some additional 
assumption such as the division of variables into classical and quantum ones are nec- 
essary (see e.g. Ref. and, as a recent contribution touching upon the problem of 



unitarity, Ref. |TB|). This type of problem has plagued early attempts to quantize 
a scalar particle in curved space-time, and similar obstructions have been encoun- 
tered when the question was raised how quantum cosmology should be interpreted. 
It might appear that no "invariant" structure of the ordinary quantum mechanical 
type — such as Hilbert spaces and unitary evolution — may be constructed at the 
fundamental level of the wave equation. This expectation has long been common 
belief, and it relies on the absence of geometric structures such as symmetries that 
would take over the role of the timelike Killing vector field in the flat Klein-Gordon 
case U (see also Ref. [T7|] for the full superspace context). 



Nevertheless, some efforts have been made to circumvent this type of problems, 
thereby making use of other than local geometric structures. Recently, proposals 
have been given that relate a preferred decomposition with the structure of super- 
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space at small scale factors |L8j and the structure of minisuperspace at large scale 
factors |n| . Furthermore, the refined algebraic quantization approach envisages the 
construction of a Hilbert space of physical states quite independent of local geo- 
metric considerations. As outlined above, it is based on the scalar product q from 
Q2.71) and thus incorporates global issues from the outset. It definitely rejects the 
interpretation of wave functions to propagate causally in minisuperspace, and hence 
the use of the Klein-Gordon type product Q from ( |2.8|) . Nevertheless, as we have 
already remarked above, its result in the flat Klein-Gordon case consists merely in 
reversing the sign of Q in the negative frequency sector. The question whether a sim- 
ilar relation exists in more general models has — to my knowledge — not yet been 
investigated systematically. All these attempts and strategies provide hints that 
in addition to the local geometry of (mini)superspace, some other structure may 
play a fundamental role in the emergence of the mathematical objects of quantum 
mechanics and the possibility of predictions. 

The aim of this article is to provide further evidence that mode decomposition 
and underlying unitarity might be achieved in a larger class of models than is usually 
believed. For the moment, our main assumption is the positivity of the potential 
U, i.e. the clear separation of modes in the classical sense. After developing the 
formalism necessary in the next two Sections, our proposed solution of the decompo- 
sition problem is formulated in Section 5. There and in Section 6 we will encounter 
hints that the construction of an underlying quantum mechanical structure along 
our proposal is related to further restrictions, possibly of the type of analyticity 
properties of the model, and that it relies on global issues. In Section 8 we will spec- 
ulate on possible relations to the refined algebraic quantization approach, and we 
give a heuristic argument how to define a preferred decomposition in this framework 
(although we have to leave open whether this decomposition agrees with the one 
following from our proposal or whether it may give rise to a different solution of the 
decomposition problem). 



3 Decompositions and unitary evolutions associ- 
ated with WKB-branches 

Our starting point is to choose a Hamilton- Jacobi function S (i.e. a solution of ( |2.5| ) 
in some domain Q) in order to define a geometric background structure with respect 
to which the wave equation may further be exploited. The question whether some 
structures are actually independent of this choice will be dealt with in Sections 4 and 
5. Hence, for the moment, we consider S to be fixed. Without loss of generality, we 
assume S a to be a future (outgoing) oriented vector. Adopting the usual notation, 
we will call S an "action" function. Since we will compute scalar products of the 
type ( f2.8| ) with £ being hypersurfaces of constant S, the domain Q shall be of the 
appropriate type. One may construct an action function by choosing any spacelike 
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hypersurface E admissible for defining the scalar product Q and imposing S = 
on E. Then the positivity of U ensures the existence of a solution S at least 
in a neighbourhood of S. The delicate problem of asymptotic spatial boundary 
conditions for wave functions will depend on the particular model. 

The function S locally generates a congruence of outgoing classical trajectories. 
The natural time parameter along these is given by the action itself. In order to 
retain the usual convention, we define the time parameter as 

t=-3, (3.1) 

future orientation corresponding to increasing values of t. (This is an immediate 
consequence of equation (|2.6|) , taking the pseudo-Riemannian signature of ds 2 and 
the positivity of the lapse function N into account). The evolution parameter t is 
not indended to be the physically experienced "time", although we take the freedom 
to use this word. The derivative along the trajctories with respect to this time 
parameter shall — in a sloppy way — be denoted as d t , and is given by 

d t = ^V a . (3.2) 

The result of the action of this derivative on a function / is sometimes written as /. 



Comparison with equation (2.6) shows that this corresponds to the choice N = U^ 1 
of the lapse function. It proves useful to introduce a coordinate system in Q adopted 
to this situation. Let £ a (a = 1, . . . (n — 1)) label the trajectories (i.e. £ a = 0) and 
choose (y a ) = (£,£ a ). It then follows that the metric is given by 

dt 2 

ds 2 = - — + lab dCde, (3.3) 

where the potential is expressed as a function U(t,£), and j a b(t, £) is the Riemannian 
(i.e. positive definite) induced metric on the surfaces Et of constant t (i.e. of constant 
action). In this coordinate system, dt is given by the partial derivative d/dt. For 
any fixed t, the functions £ a provide a coordinate system on E t . For later purpose 
let V a denote the covariant derivative with respect to r y a b. Again, for fixed t, this is 
an object intrinsic to the geometry of E t . An n-covariant version of 7 a b is provided 
by 

Paf3 = 9a? ~ = 9a/3 + ~JJ~ ■ (3.4) 

It may also be understood as the projection onto directions tangential to E t . 

We will encounter linear operators A acting on functions x(^>0 011 -M- whose 
action is actually tangential to E t (as e.g. any derivative operator composed of 
d a = d/d(, a , such as A ab (t, £)d a b)- Such an operator is characterized by 

[A, S] = , (3.5) 
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the derivative A thereof along the trajectories being defined as the commutator 
[d t , A]. (In the example just mentioned this would be A ab (t, £,)d a b). Taking again 
the commutator with d t , it follows that [A, S] = 0, i.e. A acts tangential to E 4 along 
with A. 

We are now in position to develop our formalism. In addition to S we fix a real 
function D satisfying the conservation equation 

V a {D 2 S a ) = 0. (3.6) 

In terms of the coordinates (t,^ a ), the general solution thereof reads 

D 2 (t,£)= f ® (3.7) 

where /(£) is an arbitrary function which is constant along each trajectory and which 
we assume to be chosen strictly positive everywhere in Q, and 7 is the determinant 
of the metric 7^. The function / represents the freedom to re-label the trajectories 
and defines the measure d n ~ 1 ^ /(£) (see Q3.14Q below). 

Given S and D, we express any wave function as 

ij = X De lS . (3.8) 

This resembles what is usually done in the semiclassical approach. Note however 
that we require S and D to satisfy the full Hamilton- Jacobi and conservation equa- 
tion, with all variables y involved. No division of degrees of freedom into classical 
and quantum has been assumed. Any wave function ip is associated with a func- 
tion x on the domain Q in which S and D exist. Some wave functions "fit" the 
background structure (S, D) in a natural way: When e tS varies much faster than D 
and Xi i> may be called a WKB-state. Although we will not develop a WKB type 
formalism here, we simply call the pair (S,D) a WKB-branch. (One could equally 
well call it Hamilton- Jacobi-branch; it encodes a congruence of classical trajectories 
whose tangent vector field is hypersurface orthogonal, endowed with a weight repre- 
sented by D. By using the notation admissible WBK-branch we emphasize that in a 
particular model some asymptotic condition constraining the set of possible action 
functions S might be appropriate). 

Inserting (RT5T) into (|2.2|), the wave equation becomes equivalent to 



idtX=(\d tt + h) X , (3.9) 

where h is an operator acting tangential to E t , i.e. [h, S] = 0. This may be shown 
either by direct computation or, more elegantly, by using an ADM type (n — 1) + 1- 
decomposition of the Laplacian Vo,V a . In Ref. 0] the differential equation for x 
has been worked out for more general evolution parameters. (Setting J\f = U^ 1 ^ 2 , 
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M a = in the formalism of Ref. |2D] reduces everything to the framework considered 



here). We just quote the result that h is given by 

h = H cS --D{D- 1 )" (3.10) 

where 

H °"--23u v ^ D ^-^vV v '7V V " D - (3U) 

The derivatives V a and T> a are understood in this expression to act on everything 
to the right of them (e.g. V a D acting on a function \ giving V a (Dx))- In Ref. [|(| 
it was argued that in a WKB-context ( |3.9| ) is well-approximated by the effective 
Schrodinger equation id t x ~ H cE \- {H cS was therefore called "effective Hamilto- 
nian" ) . 

Here we will proceed without approximation instead. We first compute the 
scalar product (|2.8| ) for two wave functions of the form (|3.8| ), (i.e. within the same 
WKB-branch (S,D)) for hypersurfaces £ = St. (Due to our assumptions, there 
exists at least an interval of values for t such that E t is admissible for computing 
Q). The orientation of the hypersurface element is fixed by <i£ a = <i£n a , where 
n a = —U~ 1 / 2 S a is the unit normal (n a n a = —1), and gE is the (positive definite) 
scalar hypersurface element (in the coordinate system (t, £ a ) we have n° < and 
dT, = d! 1 ^ 1 ^ v/7). Making use of (|3.2|) , the definition of n a and the fact that D is 
real, the result is easily seen to be 



jf ^ c£ VUD 2 [x* lX 2 + \{x\ (idtX*) + (^tXi)*X 2 )l ■ (3.12) 



By construction, it is independent of t. The first part of this expression amounts 
to introduce the (positive definite) scalar product for each hypersurface S t within a 
WKB-branch 

(Xi|X2> = (xi|X2> t = / dEVUD 2 xlX2. (3.13) 
In the coordinate system (t, £, a ) this reduces to 

(XiiX2> = / d n ~^ ^VUD 2 xlX2 = J ^-^1(0x1^0X2^0- (3-14) 
It thus has the Leibnitz rule property 

4(Xi|X2) = (Xi|X2> + (X1IX2) ■ (3.15) 

In the following we will make use of the formal hermitean conjugate of an operator 
A, defined by 

(XiMLX2> = (^ + Xi|X2>. (3-16) 
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In case of differential operators, it is found by formally carrying out the usual pro- 
cedures of partial integration and omitting all terms stemming from the (asymp- 
totic) boundary of E t . The conditions necessary to make them (and statements 
like A = A') rigorous are addressed to concrete models to which this formalism is 
applied. The complex conjugate A* of an operator A is defined by 

A*x = (A X T- (3-17) 

Complex conjugation satisfies (AB)* = A*B*, (Ax)* = A*x* and (A*) 1 " = (A 1 ")*. 
The operator h from ( |3.10|) is hermitean and real (h* = h* = h), and the same is 



true for the effecive Hamiltonian H eS from (|3.11|) as well as for the action S (as 
multiplication operator). The latter fact implies that the hermitan conjugate A^ 



and the complex conjugate A* of an operator A satisfying (|3.5| ) both act tangential 



to E t along with A. The dot operation A — * A = [d t ,A] commutes with both * and 



The wave equation in the form (j3.9|) is not of the Schrodinger type, but rather a 
differential equation of second oder in the classical evolution parameter t. However, 
let us try to assume that there exists an operator H such that any solution of the 
differential equation idtX — H\ automatically satisfies (|3.9|) . Differentiating this 
equation with respect to t and inserting id t x and d t tX m to (|3.9|), one finds that 
the operator C = iH — 2h + 2H + H 2 must annihilate all functions x subject to 
i&tX = Hx- We would like in addition that H acts tangential to E t , which implies 
Cx = for all x defined on E t , hence C = 0. Moreover, a short computation reveals 
that the space of all wave functions ip such that % evolves as i d t x = Hx is orthogonal 
under Q to its complex conjugate if H* = H* . The reasoning sketched in these few 
lines shall now be reversed. 

In view of the brief arguments given above, we write down the differential equa- 
tion for an operator H, 

iH = 2h-2H-H 2 , (3.18) 
together with two additional requirements, 

[H,S] = 0. (3.19) 

and 

H ] = H* . (3.20) 

A further requirement will be imposed later (see (|3.33|) below). As is easily shown 
using [h, S] = and W = h*, the two constraints ( |3.19|) -( j3.20D are propagated by 
the evolution fl3.18| ). Whenever they are satisfied for one particular value t = to, 
they are valid for all t. Hence, the system ( |3.18|) -( [3.20|) is consistent and well-posed. 



It is useful to introduce some more notation. For any solution H of the system 

let 

H +1=V = A + iB, (3.21) 
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where A ] = A* = A, B ] = B* = B and [A, S] = [B,S] = 0. The operators A and 
B provide the decomposition of H + 1 into its hermitean and anti-hermitean (and 
at the same time into its real and imaginary) part. These operators are uniquely 
determined to be 

A=1 + Uh + H*) B= % -(h^-H). (3.22) 

It is in particular the first of these operators that will play an important role later 
on. Moreover we note that, if if is a solution, the operator H = —{H* + 2) (which 
is identical to —{W + 2) on account of fl3.2U|) ) is a solution as well, the associated 
quantities being given by V = —V*, A = —A and B = B. 



We assume now that a solution H of (|3.18|) - ([3.20|) has been chosen. Recalling the 



motivation for considering this system of equations, we state the following property: 
Any solution of the equation 

id tX = H X (3.23) 
gives rise to a solution of the wave equation in the form (|3.9|), or by virtue of 



in the original form ( |2.2| ). This follows from the differential equation fl3.18|) . Let us 
denote the space of wave functions obtained in this way by 1H + . Accordingly, any 
solution of the equation 

id tX = -(H* + 2) X (3.24) 
gives rise to a solution of the wave equation as well. This is because H = —(H* + 2) 



satisfies the differential equation ( |3.18| ) along with H. The space of wave functions 
obtained in this way shall be denoted by 1H~ . We have thus defined two subspaces 
H ± of the total space M of wave functions (whereby we assume appropriate fall-off 
conditions to hold at the boundaries of S t , if necessary). 

Turning to the question how these three spaces are related to each other, we 
observe that M are complex conjugate to each other, i.e. (M^* = 1H T . In order 
to see this it is important to recall that we have chosen a WKB-branch (S, D) with 
respect to which any wave function ip is traced back to a function \- The relation 
between ip and x is always provided by (|3~8| ). Hence, if the complex conjugate 
ip* = x*De~ lS of a wave function ip shall be represented in terms of the WKB- 
branch (S,D), it must be rewritten as ip* = XcDe lS , with Xc — x* e ~ 2lS = X* e2lt 
being the associated function. (In contrast, the original expression ip* = x*De~ tS 
refers to the WKB-branch (-S, D) instead). Complex conjugation ip — > ip* of wave 
functions is thus represented as X — > Xc- It is now an easy task to show that x 
satisfies the evolution equation ( p.23[ ) if and only if Xc satisfies fl3.24| ). The number 



2 in the operator —(H* + 2) thus represents the separation between the two WKB- 
branches (±S,D) by a total relative phase factor e~ 2lS . This can be illustrated by 
a heuristic argument: In a semiclassical (WKB) context one may expect H to be 
"small" as compared to unity, and thus (|3.23| ) to describe a slowly varying function 
X- As a consequence, ( |3.8| ) is dominated by the phase factor and represents a WKB 
type state. Its complex conjugate is a WKB type state as well (with respect to the 
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WKB-branch (—S,D), i.e. with reversed trajectories), but the associated function 
Xc — the representation of tp* in the WKB-branch (S, D) — is now rapidly varying! 



Since the two spaces are complex conjugate to each other, whatever can be 
stated about M + and the evolution equation (|3.23|) has a counterpart in terms of 



IH and the evolution equation ( |3.24| ). It is thus natural to ask whether M + and 



M~ have a non-trivial intersection, and whether their sum spans the space M of 
all wave functions. The answer to these questions depends on the properties of the 
operator A. It proves useful to insert (|3 .21 ) into ( [3.18 ). The differential equation 
for H becomes equivalent to 

it + V 2 = l + 2h. (3.25) 

Separating real and imaginary parts, we find the system 

A = -{A,B} (3.26) 
B = A 2 -B 2 -l-2h, (3.27) 

where {A, B} denotes the anticommutator AB + BA. (Let us note en passant that 
yet another equivalent form of (|3.18|) is achieved by setting V = iiiu~ l . The 
resulting equation for the operator u is linear, 

u + (l + 2h)u = 0, (3.28) 

which is not a surprise since H stems from a linear problem. We will make use 
of this result later on). Equation ( [3.26D may formally be solved: Let Ao be the 
action of A on some "initial" hypersurface S to . Since A acts tangential to all E t 
(i.e. [A, S] =0), one may imagine to represent A as a t-dependent operator acting 
on functions x{0> an d then insert t = t - This makes Ao an operator for functions 
x{y)i i-e. at the same footing as all other operators considered so far. The solution 
to ( |3.26| ) is given by 

A = UA a U ] (3.29) 

where U is the operator satisfying the differential equation U = —BU and the initial 
condition IA$ = 1. Its hermitean conjugate thus satisfies W = —WB, and (|3.29|) is 



easily seen to be the unique solution to ( |3.26|) , once B is known. Hence, whenever 



there is a function xo on some S 4o such that AqXo = on E to , it follows that there 
is a function x hi the region on which A exists, such that Ax = 0. This ensures 
that there exist solutions H such that A is invertible (or positive). As an example, 
one may choose H = 0. 



There is another way to make this explicit. If ip G 1H + , i.e. x satisfies ( |3.23| ) 



then = Ax satisfies id t fi = H^fi. Likewise, if ij) G M , i.e. x satisfies ( |3.24|) 



then /i = Ax satisfies id t fi = —(H + 2)fi. If AqXo = on a hypersurface S 4o , we 
promote Xo i n t° a solution of (|3.23|) . Then fi = Ax vanishes on £ <0 and satisfies 
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a Schrodinger type equation. As a consequence \i — 0, i.e. Ax = everywhere. 
In other words, the existence of functions annihilated by A cannot be restricted 
to a single hypersurface E to . In this sense, the character of A as an invertible (or 
positive) operator is propagated by the evolution provided by ( |3.18| ). 



Subtracting ( p.24|) from ( |3.23|) , we infer that a common element of both M + and 



M~ satisfies Ax = 0. We require that this equation has no non-trivial solution. In a 
particular model, one must refer to the asymptotic boundary conditions defining M 
in order to exclude the existence of a non-trivial x- We assume that this can be done 
and end up with 1H + n M = {0}. Consequently, one may consider the direct sum 
1H + © M and ask whether it coincides with M . Given any if> G M, we consider its 
"initial data" (x,x) a ^ some hypersurface £ to . if) is the sum of two wave functions 
^ G M ± if the system of equations x + + X~ = Xi Hx + — (H* + 2)x~ = ix ma Y 
be solved with respect to x ± on some S to . The propagation of x ± °ff ^t is then 
provided by ( 3.23 ) and (|3.24|) . Eliminating x + from these equations, one arrives at 



the single equation 2Ax~ = (H — id t )x, which has to be solved with respect to 
X~ on S to . Since x an d X ma ¥ arbitrarily be prescribed on S 4o (up to the required 
fall-off behaviour at the asymptotic boundary), the decomposition of if) into a sum 
of , ± may always be achieved if A has a sufficiently well-behaved inverse. 

The structure of the scalar product for elements of depends on the properties 
of A as well. If if>i t 2 ^ JH + ) a straightforward use of ( |3.12| ) reveals 



Q{if)iM = (Xi\A X 2) , (3.30) 
whereas for ipi^ G M~ we find 

QtyiM = -bcMxz) ■ (3.3i) 

If if)± G M + and ip2 G M~ , the scalar product is 

Q(i)i,i)2) = ^(xi\(H^-H*)X2), (3.32) 



which vanishes on account of ( 3.20 ). 

In view of these results we impose an additional requirement on the solutions H 
of the system ( |3.18| )-( |3.20| ). We assume that A is a positive operator, symbolically 

A>0. (3.33) 

As a consequence, the Klein-Gordon type scalar product Q is positive (negative) 
definite on 1H + (1H~), and the total space of wave functions is a direct orthogonal 
sum 

M = H + © M~ . (3.34) 

Thus, if all steps we have presented in a rather formal way go through in a particular 
model, any choice of a WKB-branch (S, D) and any choice of H subject to the above 
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restrictions gives rise to two Hilbert spaces of wave functions. Note that for this 
construction to fix 1H ± uniquely, it is sufficient that the basic objects, S and H, 
exist in a neighbourhood of some hypersurface E to . The key role of H is thus to fix 
the initial data of ip^. At the boundary of the domain on which S and D exist, the 



relation (|3.8| ) between ip and x breaks down, thus leading to a singular behaviour 
of x an d H, while ip remains well-behaved. 

We have already mentioned that in a WKB-context the operator H can be con- 



sidered as much "smaller" than unity. In such a context, a solution of ( |3.23|) may 



be viewed to be a WKB type state "built" mainly around outgoing trajectories. 
Analogously, a solution of ( |3.24| ) may be viewed to be associated with incoming 



trajectories. Thus, the spaces M ± provide generalizations of the classical modes of 
"outgoing" and "incoming" evolution to the quantum level. In this sense, we can 
talk about a "mode" decomposition. Note however that, so far, there is consider- 
able freedom in defining such a decomposition. In Section 5 we will examine the 
perspectives for making it unique. 

The practical use we will now make of the positivity property of A is that there 
exists a positive square root A 1 ^ 2 and its inverse A~ 1 ^ 2 . So far, we have constructed 
two Schrodinger type equations ( |3.23p -( |3.24j ), but the evolution operators H and 



[H* + 2) are in general not hermitean. This may easily be cured by redefining the 
wave functions once more. If the WKB-branch (S, D) is fixed, equation ( |3.8| ) defines 
a one-to-one correspondence between if) and x- Having chosen H, we define for any 
wave function 

V = A 1 ' 2 X - (3.35) 

For ipi t 2 G M + we thus find 

QtyiM = (vM , (3.36) 

whereas for ipx^ G M~ we have 

Q{if>iM = -(m\m) ■ (3-37) 

For completeness we recall that Q(ipi,ip2) — if ipx G 1H + and tp2 G 1H~. Since 
( p. 35 ) applies for any state, the decomposition of a given wave function with respect 



to 1H reads r\ = i] + + rf 



The evolution equations (|3.23| )- (|3.24j ) may now easily be re-written in terms of 
rj. If ip G 1H + we find 

id t r] = )Cri, (3.38) 
which replaces ( |3.23|) , and for ip G M~ we have 



id t r] = -{lC* + 2)r], (3.39) 
which replaces Q3.24p . The operator /C appearing here is given by 

K = A 1 ' 2 HA' 1 ' 2 + % (A 1/2 ) -A- 1/2 . (3.40) 
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It is hermitean, 

/C f = /C, (3.41) 

and acts tangential to Sj, i.e. [/C, S] = 0. Note that in general lO 7^ fC*. The 
same properties hold for the evolution operator of (|3.39|) , in particular — (/C* + 2)t = 
— (/C* + 2). (The way how properties of /C carry over to properties of — (/C* + 2) 
may be inferred analogously to the case of the operators H and —{H* + 2) : If a 
wave function ip is represented in the WKB-branch (S,D) by the function 77, then 
the complex conjugate wave function ijj* is represented by 7] c = r]*e~ 2lS ). By using 
(EHTD , ( ggSp , the time derivative of the identity A = A l/2 A l/2 , and the fact that 



(A 1 / 2 )^ = A 1 / 2 , we can write down another expression for /C, 

K = A-\ + i -A- 1/2 ([A, B] + [A 1/2 , {A 1/2 Y])A- 1/2 , (3.42) 

which makes its hermiticity property ( |3.41[ ) explicit. In this formula, one may insert 
[A,B] = i[H,Hi\. 

Both (|3.38|) and ( |3.39|) represent (at least at the formal way we are treating 



things here) a unitary evolution with respect to the (negative of the) classical action 
S as time parameter. The coordinates £ a labeling the trajectories play the role 
of standard quantum mechanical variables. For a given wave function ip(y), the 
function r](t,£) may thus in a sense be considered as the "true" representation of 
the quantum state with respect to the background structure {S, D, H) . 



4 Variations between WKB-branches 

The WKB-branch (S, D) as well as the operator H used in the previous Section have 
not yet been specified but have been chosen arbitrarily. Here, we will examine what 
happens when a wave function is analyzed with respect to two WKB-branches. 

Let (S, D) and (S", D') be two infinitesimally close WKB-branches. They can 
be assumed to be defined in the same domain, which thus contains a family of 
hypersurfaces S t a well as a family of hypersurfaces E# . One may imagine these 
two families of hypersurfaces to be small deformations of each other. We will denote 
variations at fixed points of the manifold Ai by 5S = S' — S, 5D = D' — D 
etc. The background quantities are of course unchanged (in particular 5U = 0, 
and 5 commutes with V Q ). Finite variations can in principle be constructed from 
infinitesimal ones, as long as the initial and the final WKB-branch are not separated 
by some barrier of global type. 

The variation is characterized by 5S and 5D, which have to respect the Hamilton- 
Jacobi ( |2.5| ) as well as the conservation equation ( |3.6| ). These quantities are easily 
seen to be constrained by 

{5SY = (4.1) 
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It is convenient to write down these objects in terms of the original ("unperturbed") 
WKB-branch (S,D). In a coordinate system (y a ) = (£,£ a ), the first equation states 
that 5S(£) is an arbitrary function, and 5D(t,£) can be found by integration along 
the trajectories of S, the freedom to choose the initial value 5D(t ,£) corresponding 
to the expected freedom for D' . The variation of the derivative operator (|3.2| ) is 
given by 

m)=V = 1 fVa- (4-3) 

This is a derivative operator purely tangential to S 4 (i.e. [V, S] = 0) illustrating the 
change in the direction of the trajectories (and thus the variation of the hypersurfaces 
of constant action). 

Since the operator h from ( [3.10|) is uniquely defined in each WKB-branch, its 



variation is already fixed. Neverthless, straightforward application of 5 gives an awk- 
ward expression that may be simplified in a somewhat lengthy computation. On that 
way, one uses identities like P^SSp = 5S a , Sp5P a ^ = -5S a and [Sh, S] = -[h, 5S] 
(which follows from the fact that [h, S] = 0). After performing some operator al- 
gebra involving (anti) commutators with dt, one may get rid of several inconvenient 
terms stemming from the second contribution in ( ft .ICQ and containing derivatives 
of D. Denoting 

8= 5 R +l5S , (4.4) 

the result is 

6h=[h,S]+i(y-E)- l -[d u V -£} + [h,5S]d t . (4.5) 

Due to (|4.1|) , one may insert £ = (D~ 1 SD)' . Another expression for this variation 
is 

5h=[h,™]+±[h,6S\ + [h,8S\d t , (4.6) 

where one may replace [A, 55*] by [d t , [h, 6S] ] = [h,SS]'. Also note that [h, SS] = 
—'V + £. This form is particularly suitable for computing the variation of h and 
higher derivatives. The most beautiful and comprimed formula is 

5h=[h, 6 -^]+ 1 -{d t ,[h,5S}}. (4.7) 

In all these expressions, the appearance of dt as an operator acting directly on 
wave functions signals the variation of the trajectories. Although h' contains only 
derivatives with respect to the coordinates £ /a , it attains derivatives d t when re- 
written in terms of (t, £ a ). 
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Having made explicit how the basic quantities defining a WKB-branch transform, 
we turn to the case of H . In both WKB-branches, we have the system of equations 
d3.18[) - (p.20|) for H and ff' , respectively. In principle, both H and ff' may be 
chosen completely independent of each other, as long as they satisfy their respective 
equations (we just assume that they differ only by an infinitesimal operator 8H). 
One may of course take the variation of (|3.18|) — thereby using 5H = 5([d t , H]) = 
(SH)' + [V, if] - in order to obtain a linear inhomogeneous differential equation 
for SH that has to be satisfied anyway. (It states that, if H obeys ( |3.18| ), then if' 
obeys its respective differential equation as well). Also, one may vary fl3.19| ) and 
( |3.20p , whereby in the latter case one must take into account the subtlety that S 
does not commute with ^ . 

We will however restrict the freedom in choosing if and if' independently. We 
require that the mode decompositions induced by these two operators are identical, 
i.e. 1H + = 1H' + and M = M'~ . Any wave function can be re-written in terms 
of the two WKB-branches as ip = xDe lS = x'D'e lS . If ip is an element of both 
HJ + and IH' + , we encounter two evolution equations % d t x = if X and i d t >x' — H'x'. 
Expressing x an d Qt in terms of x' and d# and assuming 5H = H' — H to be 
infinitesimal, a simple calculation reveals that the variation of if must be of the 
form 

5H = [H,£] +i(y-£) +0, (4.8) 

where £ is given by (|4.4| ) and O is an operator satisfying Ox = for all x such 
that % d t x — if X- Taking into account that ( |3.19|) is valid in both WKB-branches 



(hence [5H,S] + [H,5S] = 0), it follows that O = [H,5S]{d t + iH). Thus, 5H is 
completely fixed, once ff has been chosen. (Supposing that O has an additional 
contribution O add , we must have [0 add ,5'] = and O add x = for all x such that 
idtX = Hx- Since O add acts tangential to E t , and the "initial value" of x on some 
S 4o is arbitrary, we infer that O add x = for all Xi hence O a,dd = 0). In other words, 
if and if' must be such that their (infinitesimal) difference is given by 

5H = [H, £}+i{V -£) + [ff , 5S] (d t + iH) . (4.9) 

Inserting V — £ = — [h, 5S], we obtain 

5H = [H, -E]+i[H-h, 5S] + [H, 5S] {d t + iH) . (4.10) 

This condition is equivalent to the equality M + = M' + . Complex conjugation 
implies the equality IH = M'~ . 

We have now a simple criterion whether the decompositions defined in terms of 
two overlapping WKB-branches are equal. Our idea is that the whole of the manifold 
M. is covered by all admissible WKB-branches of the previous type. Furthermore we 
assume that this set of admissible WKB-branches is connected in the sense that any 
two branches may be deformed into each other by a family of infinitesimal variations. 
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Whether this idea may rigorously be formulated will depend on the details of the 
particular model, its asymptotic structure and on how "admissible" is defined. Here, 
we assume that it is the case. 

As a consequence, one could try to choose the operator H within any WKB- 
branch such that for infinitesimal neighbours the condition (|4.10|) holds. This in 



turn gives rise to a unique decomposition. One may, for example, start with any 
decomposition of M and define H in each WKB-branch such that it generates the 
pre-assumed decomposition. Then ( [4.1 0|) will surely hold for close neighbours. This 
argument shows the existence of such constructions. However, in doing so, nothing 
would be gained. Our actual goal is to single out a preferred decomposition. 

At this step we must make a difference between the two possible interpretations of 
the wave equation Q2.1| ) as a quantized scalar particle and as a quantum cosmological 
model. In the first interpretation, the idea of a wave function propagating causally 
is natural if not necessary. This provides the main motivation for using the Klein- 
Gordon type scalar product Q rather than the L 2 -product q from ( |2.7| ) in quantum 
field theory on curved space fl2l |. On the other hand, no such feature is necessary 



in quantum cosmology. When trying to select a solution of the differential equation 
(13.181) , we will accept a result that lacks causality and locality properties and thus 



transcends the quantized particle picture. 



5 Perspectives for uniqueness: the iterative solu- 
tion 

In this Section we will tackle the problem whether a preferred decomposition M = 
JH + ®1H~ may be defined in a natural way. According to the formalism developed in 
the two foregoing Sections, we have to choose an operator H within each admissible 



WKB-branch (S,D), satisfying the differential equation Q3.18]) and the additional 
conditions ( |3.19| )-( |3.2(J| ). For any two infinitesimally close WKB-branches, the vari- 
ation shall satisfy the equation Q4.10Q . This guarantees that the decompositions 
associated with the two WKB-branches are identical. As already mentioned, we 
assume the global structure of the set of all admissible WKB-branches to imply the 
statement that, no matter which WKB-branch is used to define the decomposition, 
the result is alway the same. 

Thus, the next step in our program is to make an appropriate choice of the 
operator H, once (S,D) are given. This choice has to be "natural", and should 
thus not refer to additional objects (such as a prescribed decomposition). Since 
( p. 18 ) is a differential equation, one might think of fixing the "inital value" of H on 



some hypersurface £f . Two possible guesses might be H = or H = h at t = t . 
However, a simple computation reveals that the solutions of (|3.18|) defined by these 
choices cannot satisfy ( 4.10| ). The reason is that two WKB-branches differing merely 
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by a rigid translation t — > t + const must be considered different (with 5S = const). 
In other words, the value t = to does not have an invariant significance, and any 
choice of a solution H based on initial conditions will fail. Let us add that the 
attempt to use initial conditions for selecting H is in spirit related to the scalar 
particle quantization picture of our problem. The impossibility in doing so reflects 
the fact that no preferred decomposition may be singled out in a causal and local 
way. This may be viewed as a justification of the point of view that a unique mode 
decomposition is not the appropriate thing to look for in quantum field theory in 
a curved background. However, we are mainly interested in quantum cosmology, 
where new possibilities arise. 

We will now present a formal solution to the problem. It is "formal" in the sense 
that a power series emerges whose actual convergence (possibly in some regularized 
sense) may depend on the particular model. The procedure described in the following 
is thus understood as a general strategy whose application to concrete models or 
classes of models will raise a number of questions by its own. 

The starting point is to rewrite (|3.18| ) as 

H=h- -H 2 - - H . (5.1) 
2 2 y ' 

This equation may be iterated, i.e. the two terms H 2 and H may be replaced by the 
square and the time-derivative of h— \ H 2 — | H, followed by the same substitution 
whenever H appears in this new equation, and so forth. One obtains an expression 

H = h- -h 2 - -h + -h 3 + -{h,h}- -h + ... (5.2) 

in which the dependence of the rhs on H has disappeared. Only expressions involving 
h and its derivatives remain. Whenever this series converges (in a suitable topology 
on the space of operators) it defines a solution of (|5.1|) . Also, the constraints (|3.19|) - 



Q3.20|) will be satisfied on account of the properties of h and the operations \ * and 
d t . Surprisingly, this solution — if it exists — solves the decomposition problem, 
as we will now show. 



When iterating ( |5.1| ), it is useful to introduce a formal book-keeping parameter 
e whose actual value is 1, and which keeps track of the "order" of terms appearing 
in fl5~2]). By replacing h — > eh, H — > eH and d t — > ed t , the differential equation ( |5.1| ) 
becomes 

H = h--eH 2 --eH. (5.3) 
2 2 v ' 

This reflects the structure of the iteration and allows to think of e as a "small" 

expansion parameter. The rescaling of quantities by insertion of appropriate powers 

of e could have been achieved from the outset by using the evolution parameter 

t new = ef _ _ eS and the operators h new = e~ l h, H ncw = e^H instead of t, h and 

H, and thereafter removing the superscript new . As a consequence, the according 
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insertions in any of the formulae displayed in the following is unique and consistent. 
This consistency is crucial for our construction. Note that we do not introduce 
an e-dependence of h, as the second part of equation (|3.10|) might suggest. Thus, 
our particular way of inserting e's differs from a mere book-keeping of an h-type 
parameter, and it is actually simpler. It is designed such that in the procedure 
given below the operator h may be treated as an object by its own, without any 
reference to its internal structure, i.e. to its definition (|3.10|) -( |3.11|) . Nevertheless, 
the position of e in ( |5.3| ) indicates Planck scale quantities (such as m/rrip, where 
m is a mass scale in the matter sector of the theory) in a realistic model. In a 
semiclassical context one would expect H « h (in fact H m h rs H cS because the 
second term in (|3.10|) is smaller than the first one by a Planck-scale order ||20|| , a fact 
that is ignored by our way to insert e). In this sense, e may in fact be expected to 
trace "small" contributions, although quantities like h contain "small" parts as well 
which are not indicated by e. The simplicity of our book-keeping scheme precisely 
relies on this subtlety. 

We remark here for completeness that the insertion of e's along the lines described 
above changes the the wave equation Q3.9|) into id t x = (k e ®tt + h)x, whereas the 
definition ( |3.21| ) becomes eH + 1 =V = A + iB. The formula ( |3.4(Jp for the operator 



K, as well as the unitary evolution equations ( |3.3S| )- (l3.39| ) remain as they stand (if 
one rescales /C — > eJC). Equation ( |3.25|) takes the form ieP + V 2 = 1 + 2eh. The 



linear differential equation ( |3.28|) will be re- written below (equation (|5.16 )) 



Returning to (|5.3|), the iteration procedure generates a power series 

oo 

ff = £?V, (5-4) 

by which we mean that, for any given order p, a finite number (actually p + 1) 
of iterations generates the operators Ho, ■ ■ ■ 7~t p ■ Any further iteration leaves these 
operators unchanged (i.e. acts only on TC p +i,Tlp+2 ■ ■ ■)■ When the iteration proce- 



dure leads to a convergent result, it should be identical to (|5.4j) , evaluated at e = 1. 
There are two ways to formulate this situation in terms of sequences of well-defined 
operators. The first one is to consider the iterative scheme 



H n = 



1 TT o i 



H p+1 = h- - 6 H 2 p --eH p (5.5) 

for non-negative integer p. After an appropriate number of steps, the operators 7i p 
may be read off. The iterative solution to the differential equation (|5.3| ) is formally 
given by lim^oo H p . 



An alternative formulation is to insert (|5.4|) as an ansatz into (|5.3|) and separate 



orders of e. This generates the iterative scheme 

H = h 
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I p i . 

Kp+i = - o K q H v - q - - Up (5.6) 



for non-negative integer p. The corresponding solution H of ( |5.3| ) is given by the 
series (|5.4j) . The first few terms are 

Hi = --h 2 - -h 

2 2 

H 2 = lh 3 + ~{h,h}-~h, (5.7) 



and they clearly fit the structure of ( [5.21 ). In view of this scheme one could say that 
the series ( |5.4j ) is the (unique) solution to (\5.3j ) which is analytic in e. 

The convergence of the sequence ( |5.5|) implies the convergence of the series fl5.4|) • 
However, the converse does not seem to be necessarily true: it is conceivable that 
the limit of the sequence (H P )^L Q does not exist on account of effects produced 
at very large orders in e, while the operators 7i p do not "feel" these effects and 
combine into a convergent series. (A toy-example of such a situation is given by 
the scheme H = 1, H p+ i = 2eH p , where the result is H p = (2e) p , which does not 
converge as p — > oo if e = 1, but whose analogue to ( |5.4| ) yields Ti p = 0, hence 
H = 0. We also note that a modification of the scheme (|5.5|) by admitting some 
non-zero starting value Hq would change the higher order behaviour of the H p , 
while leaving the Ti p unchanged) . The actual convergence of linip^oo H p would thus 
imply these higher order effects to become successively negligible. The principle 
possibility of a situation in which they are not negligible becomes clear when one 
tries to define precisely what the dots in ( |5.2| ) mean. All of our considerations will 
concern statements given for some arbitrary but finite order in e. For this type of 
statements, both procedures are equivalent, but what we mean by "the" iterative 
solution H is — in case of doubt — the series (|5.4j ) rather than the limit of ( |5.5|) . 
Nevertheless, the iteration ( |5.5|) is more appealing as a fundamental issue because 
it does not even require the parameter e to be different from 1. In contrast, the 
procedure (|5.6j ) relies on a particular book-keeping of terms which is controlled by 
e. Thus, our hope is that these two procedures are in fact equivalent, the parameter 
e merely being a technical tool in computations that could as well be omitted. 

In the sense elaborated above, it is clear that the series (|5.4j) qualifies as a 
formal solution of the differential equation ( |5.3|) or, after setting e = 1, of ( |3.18|) . 



In case of convergence, it will be an actual solution. The same applies to the 
conditions ( |3.19| )-( |3.20| ). The really interesting question, to which we will now pay 
our attention, is whether it satisfies equation ( |4.10|) which guarantees the uniqueness 



of the decomposition. In order to investigate this question in a most economic way, 
we re-write the equations ( |4.6|) and ( |4.10|) by rescaling h, H and d t by a factor e, 
while leaving 5S and unchanged. The variation of h thus becomes 

$h =i h ^-§}+^ SS \ + e fi> SS \ d t > ( 5 - 8 ) 
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and the uniqueness condition ( |4.10| ) reads 



5H = [H, -£.]+i[H-h, 5S] + e [H, SS] {d t + iH) . (5.9) 



We also display the variation of the differential equation (|5.3|) 



5H = 6h- 1 -e{H,6H}+ t -e[ [h, 5S],H]+ l -e [H, (^) ] - U (SB); (5.10) 

where the variation of d t has been performed as shown in Section 4. We understand 
that 5h in this equation stands for the expression (|5.8| ). The variation of the formal 
iterative solution H may be computed either directly (by using explicit expressions 
like (|5.7|) ), or by using the iterative scheme (|5.6| ), or by simply iterating (|5.10|) 
(i.e. successively re-inserting ( p.3[ ) and ( |5.1U| ) for H and 5H). All these procedures 
generate the coefficients of 

oo 

5H = J2$'H P 6 P (5.11) 

to any desired order. The problem under consideration may now be treated using 
( |5.3|) and ( |5.8| )- (|5Tl0| ), without any further reference to the structure of h. As 



mentioned above, this is an immediate consequence of how we insert powers of e. 
Although we treat h as independent of e, the variation Sh involves e. This is just 
an artefact of how we group terms consistently. (We recall that the wave equation 
( |3.9|) is rescaled as id t x = (| e ^tt + h)x, which provides another explanation for the 
e-dependence of 5h). 

Once having the above formulae at hand, the argument we need is surprisingly 
simple. It is based on the observations that 

(i) inserting 5H from the uniqueness condition ( |5.9| ) into the variation ( |5.10| ) 



of the differential equation gives an identity (because (1579| ) guarantees that H' = 
H + 5H satisfies its respective differential equation along with H), that 

(ii) the converse of this, i.e. inserting the variation ( |5.11|) of the iterative solution 



into the uniqueness condition ( |5.9| ), or, likewise, inserting and re- inserting 5H from 
(|5.10| ) into ( |5.9| ), is just the equation we would like to prove, and that 



(in) the iterative structure of our problem implies that (i) and (ii) and essentially 
identical! Note that the statements (i) and (ii) hold for arbitrary order in e, due to 
the overall consistency of the way the various factors e have been inserted. 

Let SH£A denote the rhs of (|5.9| ) and 5Hfe^A the rhs of ( |5.10| ). We consider 



the equation 5H^^ = 5H^^, which still contains H and 5H. In all steps we 



understand that for H the iterative solution is inserted. Technically, this is achieved 
by inserting and re-inserting ( |5.3|) until the desired order in e is achieved. In the case 
of the remaining 5H's, we have two possibilities: Inserting SH = 5H^]fy gives the 
identity mentioned in observation (i), whereas inserting 5H = 5H ^ 1( ^ iteratively 
gives the equation of interest, mentioned in observation (ii), that we want to prove. 
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We proceed iteratively and begin with 0(e°). To this order, the equation 5H£^ = 
5H^Jih reduces to the identity [h, ^j] = [h, ^j]. The terms containing 5H are of 
order 0(e), hence to not contribute. In other words, no matter which of the two 
possibilities for inserting SH is chosen, we find that 5H(jr& = 5Hfc~{A is true to 
0(e ). 



Now we consider the equation <5-ff([TJ) = SH ^Jfy at 0(e). The remaining SH's 



enter only with their 0(e°) part, and at this order both possible choices give the 
same result, as we have shown above. Hence, we infer that SH£^ = SHfeJfy is 
true to 0(e). The following steps proceed exactly along the samelines. Using all 
previous orders, we may establish the equation SHft~^ = SH ^ ^ to hold to arbitrary 
order 0(e p ). The crucial point is that it is not necessary to specify whether we have 
to insert SH^j or 5H^ ^ for the remaining <5fTs, because the equality of both 
alternatives has been shown in the foregoing step. In other words: The uniqueness 
condition ( |5.9|) , with the iterative solution (|5.4j ) inserted, precisely yields ( |5.11[) . We 
have thus proven the 

Theorem: The formal solution (|5.4| ) and its variation ( |5.11|) satisfy the 
uniqueness condition ( |5.9| ) to any finite order in e. 



The explicit check thereof for the first few orders is in principle straightforward, 
but rather time-consuming. Before the general pattern of the proof has become 
clear, we have explicitly done the computation for 0(e) by hand and for 0(e 2 ) and 
0(e 3 ) using Mathematica 2.2. 

This is our proposal for solving the decomposition problem. Applying the scheme 
to concrete models, one will have to make sure that the the series fl5.4Q makes sense. 
In Section 7, a number of examples is discussed. 

We are now in position to give the first few terms in the expansion of the quan- 
tities A, B and K, which have proven to determine the unitary (Schrodinger type) 
evolution within the spaces with respect to a WKB-branch. Using the her- 
miticity property of h and its time-derivative, it is easy to compute H*. Then ( |3.22| ) 
immediately yields (taking into account the rescaled definition eH + 1 = A + iB) 

A = 1 + eh- ^e 2 /i 2 + e 3 Q/i 3 - i^+0(e 4 ) (5.12) 

B = -^e 2 h + ^e 3 {h,h} + 0(e 4 ). (5.13) 

The first of these equations could provide the key for the question whether the pos- 
itivity property (|3.33|) of A that was necessary to write down the unitary evolution 



equations ( |3.38|) - (|3.39| ) actually holds in concrete models. Proceeding formally, the 
operator A 1 ^ 2 defining the effective wave function by (|3.35|) becomes 



A^ = l + \eh- le 2 h 2 + e 3 (^h 3 - ^)+0(e 4 ), (5.14) 
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the inverse of which as well as the commutator with <9 t being readily computed. 
Inserting the results into into (|3.40 ), the evolution operator K, takes the form 



K = h- ±eh 2 + e 2 (±h 3 - l -[h,h}- i/i)+0(e 3 ), (5.15) 

whose hermiticity is evident. Whenever the position of e characterizes quantities 
that are actually small (e.g. in a semiclassical context), these expressions may be 
used as approximations. 

Let us at the end of this Section comment on an alternative formulation specify- 
ing our proposed solution (|5.4|). Recall that the differential equation for H is related 
to the linear problem (|3.28| ) which now reads 



e 2 u + (1 + 2eh)u = 0. (5.16) 

Any solution thereof povides a solution to ( |5.3| ) by V = 1 + eH = ieim~ x . It is easy 
to show that any solution u of the form 



u = v e 



~ it/e with v = ^ e<? ( 5 - 17 ) 

g=0 



provides the iterative solution Q5.4p by 

H = ivv~ 1 . (5.18) 

Here it is of course understood that the coefficient operators v q do not depend on e. 
Although the solution of the form ( |5.17| ) is not unique, all possible freedom cancels 
in (|5.18|) so as to make H unique. (If h were a c-number function of t, one could 



state that (|5.17|) is unique up to multiplication by an analyic function of e). At 



the formal level (|5.18|) reproduces the coefficients of the power series ( |5.4|) . The 



differential equation satisfied by v turns out to be 

i v — hv H — e v , (5.19) 

which may be used to compute the coefficients v q by iteratively solving the linear 
inhomogeneous differential equations iv = hv , iv q+ i = hv q+ i + (and choosing 
arbitrary integration constants). The formal series emerging from this procedure is 
nothing but the WKB-expansion associated with the differential equation ( |5.16|) to 
all orders in e (which serves as the "small" fr-type parameter). This analogy becomes 



even more appealing by the fact that (|5.16|) formally resembles a one-dimensional 



time-independent Schrodinger equation, with t being the position coodinate and 
— 1 — 2eh playing the role of the potential. In a sense, the computation of H is 
traced back to a problem within the framework of conventional quantum mechanics. 
We may thus be confident that in a large number of concrete models our proposed 
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solution H actually exists (i.e. that the formal power series for v and H belong to 
analytic operator valued functions). 

In the literature, the general solution of a differential equation of the type u + 
Vu = is sometimes given as a linear combination ve~ lt + we lt or r sint + scost. 
However, the requirement that u is of the form (|5.17f) is not an attempt to say 



something about the time-dependence of u for finite e, neither is it an asymptotic 
condition for large \t\ mimicking "positive frequency" (we do not even suppose the 
coordinate t to extend over all real values). There is no fall-off condition for v 
involved. It may happen, for example, that for some finite e the oscillatory factor 
e -*t/e j n u j g num erically almost irrelevant as compared to v. An illustration of this 
type of solution is provided by Example 2 in Section 7 where H is computed for 
h = a + bt and v turns out to behave oscillatory and bounded in the domain h ^> 
but oscillatory and exponentially unbounded in the domain h <C 0. The nature 



of v as defined in ( 5.17 ) without any further condition is absolutely crucial in our 
construction: for generic time-dependence of h, one would not expect a solution of 
the form ( |5.171 ) to exist if v was supposed to be "small" or "slowly varying". The 
behaviour of the quantities u, v and H should not be confused with that of wave 
function (despite the fact that ( |5.19| ) is of the same form as the wave equation (|3.9|) ). 



There is another interesting feature related to ( |5.17|) . Suppose the time depen- 
dence of h is C^°, i.e. h vanishing identically outside an interval to < t < ti but 



being non-zero inside and C°° for all t. As a consequence, u from ( |5.17| ) is purely 



"positive frequency" outside the interval (u = c±e %t l e for t < t and t > t\ with 
c± being independent of t). On the other hand, ( |5.17| ) is a differential equation 



that can be integated once u = c_e~ rf / e for t < to is known, and the result will in 
general be such that u ^ c + e~~ %t ^ for t > t\. This tells us that there will be no 
solution u of the required type. (This is actually the argument against the existence 
of a preferred standard-type separation into positive and negative frequencies in the 
non-stationary case). Since v and H may in any case be generated as formal power 
series, we conclude that in the case of h being in t, these formal series will not 
converge and cannot be associated with analytic functions in a unique way. Similar 
effects can be expected when the basic ingredients of the model (ds 2 and U) involve 
quantities. It is not quite clear what the correct conditions on the model are 
in order to make sense of the series (|5.4j) , but it is conceivable that they have to do 
with analyticity. 



6 Discussion & Speculations 

Let us make some speculations concerning the structure we have possibly touched 
upon. One appealing feature of the iteration procedure defining H (preferably in 
the version ( |5.5| ), although version ( |5.6| ) is technically simpler) is that a (formal) 
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solution is specified "without ever having made a choice". The procedure as such 
just consists in shifting the unwanted if -dependence at the rhs of ( |5.3|) to arbitrarily 
large orders in e. The key criterion for solving the problem is the uniqueness condi- 
tion ( |5.9| ). Nevertheless, one would like to know what is behind such a procedure. 
The whole situation is a bit remniscent of the idea of general covariance, stating that 
the geometric objects used to formulate physical laws can be described in terms of 
coordinates, but still have their right own as "objects". The different descriptions 
(e.g. of a tensor field in various coordinates) must match (i.e. obey appropriate 
transformation laws) so as to allow for an "object" -interpretation (e.g. as a ten- 
sor). Considering this as an analogy to our situation, one is tempted to talk about 
"covariance with respect to WKB-branches" . The WKB-branch {S,D) thus plays 
the role of a coordinate description of a geometric structure. The corresponding 
structure in our case encodes not only the (unique) decomposition, but also leads to 
the unitary evolution equations ( 3.38 ) and ( 13.391 ) within the two subspaces IH ± of 



wave functions. These equations show up only in the context of WKB-branches, i.e. 
there is no unique hermitean operator, just as a scalar field is not given by a unique 
function of its coordinates unless the coordinate system is specified. Nevertheless, 
the structures emerging in WKB-branches, such as decomposition and unitarity, 
match in a beautiful way. 

We should emphasize here that we have not dealt with the question how the 
evolution parameters t relate to what we experience as time, nor have we sug- 
gested a complete scheme how to reconstruct standard quantum physics and "the" 
Schrodinger equation. We merely suggest that, whenever the decomposition problem 
may be solved positively (i.e. the series for H can be given a sense) the mathematical 
structures of standard quantum mechanics (positive definite scalar product, Hilbert 
space, unitarity and the possibility of writing down probability-type expressions) 
begin to come into our reach. It is thus conceivable that the further steps necessary 
in deriving quantum mechanics from the Wheeler-DeWitt equation proceed within 
this familiar mathematical framework. 

The most difficult of the questions raised seem to be the technical aspects of 
the convergence properties of the formal solution ( |5.4|) . Without trying to go into 
the details here, we note that time-derivatives of h of arbitrarily large order are 
involved. Thus we expect the solution H, whenever it exists, to rely on the global 
structure of the model and the WKB-branch. As we have pointed out at the end 
of the preceding Section it might be some kind of analyticity condition that ensures 
convergence of the series. We have seen that at least C°° is likely not to be the 
appropriate differentiability class of (ds 2 ,U). On the other hand, it is conceivable 
that the need for analyticity is just a technical tool at some stage in order to give 
the operator H as provided by the formal series ( |5.4| ) a definite meaning, and that 
this condition disappears when the action of H is maximally extended. This might 
be in analogy to the attempt to define the operator hi = exp(e^) by means of its 
formal power series. As an intermediate step, one encounters analytic functions, but 
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in the end, the action of U may be extended to the whole of the Hilbert space L 2 (M), 
without any reference to analyticity. If this analogy holds, the structure responsible 
for the preference of a decomposition would be of global (geometric) nature, rather 
than relying on analyticity, and the class of models to which our formalism applies 
can be expected to be rather large. 

In any case, our proposal seems to transcend the standard local differential geo- 
metric framework (within which we know the solution of the generic decomposition 
problem to be negative 13). As we will illustrate for the case of a toy-model in the 



next Section (Example 7), a naive treatment of the series (5.2) or ( |5.4|) leads to an 
argument against the existence of the unique decomposition that might be erraneous 
though intuitively appealing. Maybe one could state that, despite the lack of a local 
geometric structure (such as a symmetry) specifying a unique decomposition, some 
underlying global or analytic structure does the job. In the latter case our proposal 
might give a hint towards a subtle role, not yet completely understood, played by 
analyticity in the proper formulation of the laws of nature at the fundamental level. 



The major difference between the picture suggested by a quantized scalar particle 
and the picture due to a quantum cosmological model has already been addressed 
above: In the former, one thinks about the wave function ip as of a field propagating 
causally in space-time, and the interpretational framework, such as a decomposition 
of the space of wave functions, should be in accoordance with this point of view. 
Thus, whatever the significance of such a decomposition is, it should be determined 
locally on (or near) some spacelike hypersurface and not depend on what will hap- 
pen to the background (Ai,ds 2 , U) in the far future. The local character of the 
standard frequency decomposition in generic backgrounds refers to such a situation. 
In contrast, the appearance of a formal series in our proposal destroys this causal 
picture. The decomposition of M into M ± is determined by global issues in that it 
cannot be inferred from knowing (ds 2 , U) only near some hypersurface. The scalar 
particle picture can at best lead to locally defined decompositions, e.g. associated 
with spacelike hypersurfaces. It is conceivable that our proposed unique decompo- 
sition is related to all these local decompositions by some kind of average. (We 
can even imagine that the actual well-posedness of such an average is tied to the 
convergence of the series for H. This would trace back the convergence problem to 
the existence of, say, an appropriate measure). Since the idea of causal propagation 
of the wave function is no longer crucial in quantum cosmology, these considera- 
tions are not necessarily objections against our approach but might help to uncover 
the structure which it is based upon. Also, there might be relations to the refined 
algebraic quantization scheme, on which we will speculate in Section 8. There we 
will give a heuristic argument that this approach seems to require the existence of a 
preferred decomposition (whose relation to ours remains open). 

In order to get some feeling for the way how the formal solution may become an 
actual one, we consider the simple case in which h = 0. In terms of the coordinate 
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system (t, £ a ) this means that the operator h does not depend on t. The flat Klein- 
Gordon equation (ds 2 = —dr 2 + dx 2 , U = m 2 , S = —mr = —t, D — 1) is a special 
case, with h = -^jA. (We note that D = const is not a general consequence of 
h = 0. In two-dimensional models, n — 2, one can show that h = implies D — 0, 
but even in this case .D may depend on £). For a general situation with h = 0, the 
iteration procedure generates a formal solution with H = 0. Thus, ( |5.3| ) reduces to 
a quadratic equation for H. One may easily check that the series (|5.4|) 



H = h--eh 2 +-e 2 h i - -e 3 h 4 + 0(e 4 ) (6.1) 
2 2 8 



is just the power series of the closed expression 



H = - [VT+2eh - 1) (6.2) 



around e = 0. Since one has to insert e = 1, this expression makes sense if h — \ is a 
non-negative operator, the symbol being understood as providing the (unique) 
non- negative square root. In the case of the flat Klein-Gordon equation ( |6.2|) leads 
to the expression mH = [m 2 — A) 1 / 2 — m, which was already mentioned in Section 2 
(and denoted E there). It may be well-defined in terms of the Fourier representation, 
where —A just becomes k 2 . 

In more general models one will encounter generalizations of this positivity con- 
dition. We do not expect these to create major problems, and as a hint we note that 
the condition h — | > is likely to hold in realistic models. (In the semiclassical con- 
text one even has h 2 1; see Ref. PU| ). In general, the structure of ( |3 . 1 0| ) — ( |3 . 1 1| ) 



tells us that H eS is (at least at the formal level at which we are treating such issues 
here) a non-negative operator, and whenever the "potential term" — ^D^D^ 1 )" in 
?.10|) is bounded from below, we can expect h to share this property. 



The results obtained for the case h = shed new light on the nature of the 
parameter e and the way how to re- insert e = 1. If h were a number, the series 
( |6.1|) would converge only if \2eh\ < 1. However, the symbol h actually stands for 
all numbers contained in the spectrum of the operator h. This situation may be 
illustrated in the flat Klein- Gordon case h = — ^^A. Applying h on a function 

X (x)= Jdke**x$), (6.3) 

it effectively acts as i^k 2 . Thus, the closed expression ( |6.2j ) leads to the appearance 
of \Jl + ek 2 /m 2 . The power series thereof converges only if \ek 2 \ < m 2 . Thus, we 
may apply the series ( |6.1| ) to a function x(%) whose Fourier transform x{k) has 
support inside a bounded interval — k < k < k . Choosing < e < m 2 /fcg, the 
series converges and argees with (|6.1| ). Although e = 1 may lie outside the domain 
of convergence, it is always possible to analytically continue the result unambigously 
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along the real line to e = 1. This procedure may be performed for all functions with 
bounded support in momentum space. Since the set of these functions is dense in 
any reasonable topology considered for this situation, the action of (|6.2| ) may be 



reconstructed from the series ( |6.1| ) although it will not converge on any function. 
Thus, it might be the case that in more general situations the structure of the series 
( |5.4|) should first be analyzed for small e, restricted to an appropriate set of states 
for which it converges. The original value e = 1 would then be achieved by means 
of analytic continuation. 

The above considerations seem to impy a suggestion how to proceed in more 
general cases. One could define, as a preliminary space, the set of all wave functions 
X for which there exists a non-zero e (and hence a non-trivial interval of e's) such 
that 

oo 

E^(^X) (6.4) 

p=0 

converges in a reasonable sense. The action of the power series operator (|5.4|) 
should then be constructed by analytic continuation and an appropriate extension 
to (diskrete or continuous) linear combinations. This is of course only a schematic 
sketch, and in practice various technical problems may occur. 

In case the series ( |6.4|) for H turns out not to converge on a sufficiently large 
number of states (or even on none), one might think of applying a regularization 
scheme. In this case one would have to check the branch-independence of the reg- 
ularized decomposition. One might also speculate on the existence of some mathe- 
matical structure that is poorly accounted for by our power series for H, with the 
effect that the first few terms give reasonable results, whereas the thing as a whole 
diverges. (Note that a similar issue arises for the perturbation expansion in ordinary 
quantum field theory). A minimalist point of view would, by the way, consider the 
operator K, as the most important object for practical computations. In case there 



are models in which its formal power series (i.e. (|5.15|) to all orders) makes sense 



whereas (|5.4j ) does not, one might be tempted to accept non-regularizable infinities 
in H. 

Let us close this Section by another speculative remark. We can imagine that 
procedures other than the one leading to ( |5.4[ ) exist and solve the decomposition 
problem at a formal level. (Such procedured could e.g. be based on differently 
defined iteration schemes, due to alternative way to group terms by means of a 
book-keeping parameter e) . It is not clear whether different schemes effectively lead 
to the same result. In case they do not, all of these non-equivalent schemes would 
provide a setting within which the "natural" construction of a unique decomposition 
is possible. They would be regarded as different realizations of what was called "co- 
variance with respect to WKB-branches" above (just as different tensor fields might 
solve some geometric theory), and should be considered as inequivalent solutions to 
the decomposition problem. The naturalness of the construction of a decomposition 
(lacking reference to external structures) may be regarded to bear some analogy to, 
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say, the problem of finding tensor fields built merely from a metric. Since the actual 
existence of our proposed solution seems to rely on further conditions, as discussed 
above, the number of different "natural" solutions might depend on the model as 
well. Whether all these may in principle be "realized" (just as different metrics 
solving Einstein's field equations may realize space-time) or whether one of them is 
singled out fundamentally, is of course at the moment impossible to answer for us. 



7 Examples 

Whenever a WKB-branch (S, D) has been chosen and the operator h is known, 
the key procedure in our framework is to evaluate the formal solution fl5.4|) of the 
differential equation ( |5.3| ) and to determine whether (and on which functions) it 
actually converges. This justifies the consideration of ( |5.3j ) in situations in which 
the operator h is simple. The property of h being a second order differential operator 
in the coordinates £ a may even be ignored to some extent, and emphasis may be 
laid on its dependence on the evolution parameter t. In the following we list some 
simple cases that can be treated exactly, and a scenario illustrating a class of models 
within which our proposal may directly be confronted with the standard frequency 
decomposition picture. 

Example 1: h — 

This it the simplest case, and it has already been studied in Section 6. The 
power series ( |5.4|) corresponds to the closed expression ( |6.2|) . Due to the simplicity 
of this situation, non-commutativity of operators does not play any role here. In 
the case of the flat Klein-Gordon equation h is proportional to —A, and H may be 
given a definite meaning, as outlined in Section 6. 

Example 2: h — and [h, h] = 

This is the special case h = a + bt, where the operators a and b do not depend 
on t and commute with each other. We will treat them as numbers, and in par- 
ticular assume that b is invertible and \b\ exists. The general solution of the linear 
differential equation (|5.16|) can be expressed in terms of Airy functions (Ref. 



p. 446). Using their asymptotic expansion formulae, the particular solution of the 
form ( |5.17| ) turns out to be (up to an e-dependent factor which is irrelevant) 

u = Bi(— x) — i sgn(/i) Ai(— x) (7.1) 

with 

1 + 2eh , . 

x = — —. -. 7.2 

(2e 2 |/i|) 2 /3 V ; 

Hence, the proposed unique solution for H is given by 

ie du (l + 2eh) 1 / 2 du .„ . 
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This results into the corresponding expansion ( |5.4j ) in terms of e, the various sgn(h) 
and \h\ terms being properly absorbed into an expression analytic in h (which thus 
admits h = as a special case). An alternative way to expand H is to keep the 
combination 1 + 2eh and treat only the e in the denominator of ( |7.2|) as expansion 
parameter. (This actually amounts to re-interpret the differential equation ( |5.3|) by 
distinguishing the two e's as two different parameters and expanding with respect 
to the second one). This type of expansion is given by 

eH = (VT+2^h - l) - - + - . ^ h) " + 0((e 2 h) 3 ) , (7.4) 

V > 2 l + 2e/i 8 (l + 2e/i) 5 /2 vv ; ; ' v ' 

thus providing a generalization of ( |6.2| ). The factor — | of the term linear in h may be 
checked against the second term of Tii in the general expansion (|5.71 ) . The existence 
of this solution shows that the series (|5.4|), computed by means of the iteration ( |5.6| ) 
with h and all higher derivatives neglected and operator ordering being ignored, is 
actually the power expansion of an analytic function. 

This example provides an illustration that the form ( p. 17\ ) is not indented to 
impose conditions on the asymptotic behaviour of v, because u blows up exponen- 
tially in the domain ft<0 while it is bounded and oscillates for h ^> 0. Hence, the 
function v is oscillating but has an exponential prefactor in the domain ht < 0. 

Example 3: h = polynomial in t with commuting coefficients 

This case seems to be untractable in general, but in principle works like the 
previous one when again use of the linear differential equation ( |5.16| ) is made. The 



case of h being quadratic in t leads to parabolic cylinder functions (Ref. f22f 
685). In general, the power series of the proposed solution H coincides with the 
result of the iteration procedure ( |5.6|) when all derivatives of h higher than the 
degree of the polynimial P are dropped. 

Example 4: h = a 

This case is distinguished by its scaling properties: a glance at the series (|5. 1| ) 
shows that H is e times a function of the ratio t/e. The linear differential equation 
( |5.16|) leads to a special case of confluent hypergeometric functions, the solution of 



the form ( |5.17| ) being (up to an irrelevant non-integer power of e) 



u = U(ia, 0, z) e~ lt/e with z = ^t. (7.5) 

Here, U is the second Kummer function (Ref. f2~2fl , p. 503) 

00 (a) (a + 1) (—) q 
U{a,0,z) =z- a F{a,z) with F{a,z) = J2 1 M d£ L_L ; (7.6) 

where (6)o = 1, (b) q = b(b+l) . . . (b+q—1) for any b. The proposed solution of ( |5.3j ) is 
given by H = iiiu^—e^ 1 = ivv^ 1 with v — as defined in Q5.17D — being proportional 
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to U(ia, 0, z). Using the form of u as displayed above, we get H = ~ — 2 -F~ 1 d z F 
with z = — t inserted. Thus, H is analytic in e and in fact coincides with the series 
(|5.14|) for h = at~ x . Denoting its functional dependence by H(a,e,t), we have the 



symmetries H(a,—e,—t) = —H(a,e,t) and H(a,—e,t) = —H*(—a,e,t) (all three 
variables are assumed to be real), hence H(—a, e, —t) = H*(a, e, t). For large \t\ the 
ratio H/h approaches unity. 

Example 5: h = f3t~ 2 

This case is particularly important because it emerges in simple quantum cos- 
mological minisuperspace models. In the spatially flat Friedmann-Robertson- Walker 
(FRW) universe with scale factor a (being the only minisuperspace variable) and 
positive cosmological constant A, the DeWitt metric and the potential are given by 
ds 2 = —da 2 and U = a 4 A. The WKB-branch is specified by S = — |a 3 V / A = —t and 
D = t -1 / 3 (it is unique up to S — > ±S + const and D — > const D), the range of the 
"time" coordinate being t > 0. Equation (|3.10|) gives (setting H eS = 0) h = |r 2 . 

A slightly more sophisticated model is obtained by putting a massless scalar 
field into the FRW universe, thus ds 2 = —da 2 + a 2 dcf) 2 and U = a 4 A. The 
action is chosen as before S = — |a 3 \/A = —t, whereas D = t^ 1 ^ 2 . The variable 



<p plays the role of the coordinates £ a , and applying ( |3.10| )- (|3.11|) one easily finds 
h — 1 1~ 2 (— | + |). For any t > this is a positive operator on the hypersurface 
Ti t (i.e. acting on functions x(0))- m more realistic models it is often extremely 
difficult to specify an exact solution of the Hamilton- Jacobi equation, but we can 
expect a t~ 2 behaviour in h to play a dominant role as well. 



The general solution of the linear problem ( 5.16j ) for h = (5 1 2 may be expressed 



in terms of Bessel functions (Ref. ||22| , p. 358). Confining ourselves to the range 
t > for the moment, the combination involving the second Hankel function 

u = ^Hi 2 \z) = ^(J u { z )-tY v {z)) (7.7) 

with 

z — - and v — 7.8 

e 4 e v ' 

is of the form ( |5.17f) , up to an irrevevant e-dependent factor. Expanding u for small 
positive e at constant v and positive t (i.e. using the asymptotic expansion of Bessel 
functions for z — > oo) gives an expression of the structure cve~ %t ^ where c depends 
on v and e but not on t (and is thus irrelevant for our purposes) and v is a power 
series in e that depends on v only through v 2 . Inserting the second equation of 
( |7.8| ) into v yields a power series in e. Note that j3 may be of either sign or zero, 
and v will be real or imaginary, depending on the values of (3 and e. The proposed 



solution for H is given by if = iiiu^ 1 — e _1 or, likewise, by ( 5.18 ), its expansion 
in powers of e reproducing the formal solution ( |5.4|) for h = (3t~ 2 . In the closed 
expression for H, any of the two solutions for v may be inserted because the ratio 
H^\z)~ l d z H^\z) is invariant under the substitution v — > —v. For fixed positive 
e, the ratio H/h approaches unity in the limit t — > oo. Due to singluarities in 
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the Bessel functions, the closed expression is convenient only for positive t and e. 
Nevertheless, the power series reveals that H is analytic in j3, e and t~ l . Denoting 
this dependence by H((3,e,t), the other ranges of variables are provided by the 
identities H(f3, e, —t) = H*(f3, e, t) and H({3, -e, t) = —H*(—f3, e, t) where /3, e and t 
are understood to be real. The simultaneous change of sign of these three variables 
changes the sign of H. 

Example 6: h = k + (3t~ 2 

This case may be treated in complete analogy to the previous one. The relevant 
solution of Q5.16Q is given by 



u = VtHj y 2 \ct) with c = -Vl + 2ek (7.9) 

and v as before. The preferred decomposition is therefore defined by H = iiiu' 1 — 
e _1 . The factor c replacing e -1 in the argument of the Hankel function has the effect 
already encounted in the first example that the ratio H/h does not aproach unity 
but (e£;)- 1 (v / l + 2e£;-l) for large t. 

A further possible generalization is h — k + at~ 1 +/3t~ 2 which leads to Whittaker 
functions (Ref. pj, p. 505). 



Example 7: h being asymptotically constant in two regions 

Here we would like to discuss an argument that might be raised against the exis- 
tence of a preferred decomposition. Suppose h is a c-number function of t converging 
to a constant as \t\ — > oo but undergoing some non-trivial smooth change at fi- 
nite times. For simplicity we assume hoo = and set D = 1. Naively, one would 
expect the form of the series (|5.4j) to indicate that H vanishes asymptotivally for 
large \t\. If this were true, the proposed decomposition could be thought of agree- 
ing with the standard asymptotic flat space-type decomposition into positive and 
negative frequencies {ip ~ e Tlt ) in the "far past" region t — > — oo and at the same 
time it would agree with the according frequency decomposition (again ^ ~ e Tlt ) 
in the "far future" region t — > oo. However, we know that in case of a nontrivial 
time-dependence of h in some intermediate region, these two frequency decompo- 
sitions are different. This is because a wave function starting with, say, positive 
frequency in the far past (ip ~ e~ lt ) will "feel" the time-dependence of h and in 
general evolve into a linear combination of both frequencies e Tlt at late times. (In 
terms of the physics of a scalar field, this indicates the production of particles by 
the time-dependent background. In terms of mathematics, it is represented by a 
Bogoljubov transformation). Thus, the argument states that the proposed preferred 
decomposition coincides with the two asymptotic ones which are however different: 
this is obviously a contradiction. The answer is presumably that the expectation 
H — > in the asymptotic regions is in general not justified. In the Examples 4, 
5 and 6 treated above, h vanishes for large \t\ but has a singularity at t = 0, so 
that these cases are no good illustrations for the issue under consideration here. 
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A better example would be h = /3(1 + t 2 /a 2 )~ l . In this case either H does not 
exist (its formal series (|5.4j ) diverging for all e 7^ and for all t) or does not vanish 
asymptotically. In any case, the non-local character of the principle specifying H 
is clear. Although having no explicit example at hand, these considerations should 
illustrate how our proposal might be capable to transcent the point of view that a 
preferred decomposition necessarily relies on local geometric structures. In view of 
the speculations made in Section 6, the preferred decomposition — in case it exists 
- might turn out to be an average over the asymptotic frequency decompositions. 

8 Relation to refined algebraic quantization? 

The results achieved in this paper have been of a rather formal nature, and the pre- 
cise condition under which they materialize into well-defined mathematical objects 
is not known (apart from the speculation that it refers to the global structure of 
minisuperspace and possibly has something to do with analyticity). This makes it 
difficult to re-express our findings in terms of an underlying structure or object for 
which the equations we wrote down are just representations with respect to local 
WKB-branches. 

Nevertheless, in pursuing the subject, one might gain some unified view of what 
happens when different quantization and interpretation schemes are applied. As 
already mentioned in Section 2, the interpretational framework associated with the 
wave equation (|2.2|) can be based either on the L 2 -scalar product q from ( |2~7l) or 
on the Klein-Gordon type scalar product Q from ( |2.8| ), the latter being particularly 
appropriate for the scalar particle quantization picture. Our approach is in some 
sense a mixture between the use of Q (which admits the notion of causal propagation 
on the manifold M) and a proposal for solving the decomposition problem that is 
motivated by quantum cosmology and transcends the notion of causality. 

This opens perspectives for a deeper understanding of the relation between vari- 
ous structures associated with the wave equation (|2~2j ). With regards to our proposal 
for fixing a decomposition M = M + © H~ , the most natural question seems to be 
whether the space (M,Q phyB ) = (H + ,Q) © (]H~,—Q) is closely related (or even 
identical to) the Hilbert space of physical states (TC ph s , ( | ) phys ) as constructed by 
the refined algebraic quantization approach, i.e. to what extent the results achieved 
for the flat Klein-Gordon equation [/J, as mentioned in Section 2, carry over 
to more general models. In the positive case, our approach seems to provide a piece 
of information that does not show up explicitly in the refined algebraic scheme, 
namely the decomposition, i.e. the natural appearance of two orthogonal subspaces 
M ± spanning the whole of M. Even if the answer to this question is more compli- 
cated, one can reasonably expect to learn more about the relation between Q, q and 
( I ) phys . Eventually, one might be able to reconcile the different approaches in terms 
of a single, unified understanding of quantization. 
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Although we have not entered into the formalism of the refined algebraic quan- 
tization, there is a heuristic argument relating the inner (positive definite) product 
( I ) P hys to the (indefinite) Klein-Gordon type scalar product Q. (A construction of 
the type we will now describe seems to have been written down first in the contex 



of de Sitter space by Nachtmann long ago |TTJ . In a more general context a similar 
reasoning appears in the work of Rumpf und Urbantke [I2[). Although the Hilbert 
space 7i h e of physical states emerging in this quantization scheme is constructed 



in a rather abstract way, is seems likely P3 that it can be represented as a set of 



solutions to the wave equation ( |2.2|) — modulo the issue of completion — , i.e. that 
it essentially agrees with M. The space M thus carries two scalar products, ( | ) phys 
and Q. We may find some basis {S r } of wave functions which are othonormal with 
respect to the former scalar product, i.e. (S r |S s ) phys = 5 rs . Generic states may be 
expanded as ip = J2r ipr^r- The matrix K rs = Q(E r , S s ) defines a linear operator K 
by (Kijj) r = J2s Krsips which relates the two scalar products by 

QfoMHVWW- (8-1) 

The properties of Q imply that the matrix K rs (and hence the operator K) is 
hermitean with respect to ( | ) phya and invertible. (Actually K is invertible if Q is 
non-degenerate, i.e. if Q(ip,<f)) = for all ip implies = 0. There are situations 
(ds 2 , U) in which this is not true, i.e. in which harmless initial data (^|s, ^"VaV'ls) 
may lead to an exponential blow-up with time and thus have to be excluded from 



M. Accordingly, Q(?/>, </>) = for all ip G ffl with 0^0 becomes possible |T] . 
Athough this is not likely to happen in reasonable quantum cosmological models, it 
restricts the number of mathematically admissible models for which our argument 
applies. Hence, we assume the model under consideration to be such that Q is non- 
degenerate). If an operator K defined by (|8.1|), it is unique because (i/;\K(f>) , is 
known for all ip and 0, and its matrix elements with respect to the basis {S r } are just 
K rs . Since K is invertible, its spectrum decays into a positive and a negative part. 
(Thereby we actually assume that K is self-adjoint with respect to ( | ) . , a point 
that should be made more rigorous in a deeper investigation). The space M may 
thus be expected to decompose into two subspaces IH ± , associated with the positive 
and negative part of the spectrum, respectively. In other words, is spanned by 
the (generalized) eigenvectors of K to positive /negative (generalized) eigenvalues 
(the latter appearing in pairs (A, —A)). In addition, there are good chances for 
the spectrum of K to be bounded off zero. (For the flat Klein-Gordon equation 
only one pair (A, — A) with A ^ occurs). The two supspaces M ± are orthogonal 
to each other with respect to ( | ) phys and hence, using (^.1[), also with respect to 
Q. Moreover, it is clear that if this construction can be made rigorous, complex 
conjugation maps M + onto M~ and vice versa in a bijective way. This argument 
shows that the inner product ( | ) phys is likely to serve as a structure in addition to 
the local geometry of minisuperspace and to single out a preferred decomposition. 
Since the procedure leading to ( | ) phya in the refined algebraic quantization scheme 



is essentially unique, the decomposition derived from (|8.1| ) should be unique as well 



39 



Hence, it should rely on the "geometry" of minisuperspace in some global sense. It 
is however unclear whether this decomposition coincides with the one proposed in 
Section 5, i.e. the one based on the formal solution ( |5.4|) for H. 

This point — at least as a theoretical possibility — seems to have escaped the 
attention of several authors that have contributed to the development of the re- 
fined algebraic quantization approach. For example, Higuchi || quotes Kuchaf's 
discussion ||, |17| that the absence of symmetries of the Wheeler-DeWitt equation 



does not allow a unique separation of positive- and negative-frequency solutions 
and simply states that his method (which is essentially the refined algebraic ap- 
proach) "does not require such a separation" . Also, Ashtekar et al || remark that 
the construction ofTi, (in the context of the Klein-Gordon equation on a curved 

pnys \ - 1 - 

background space-time) "may come as a surprise to some readers, as it seems to 
violate the accepted idea that there is no well-defined notion of a single particle in a 
non- stationary space-time ". Their resolution of this paradox is that the new quan- 
tum theory "contains no notion of a conserved probability associated with Cauchy 
surfaces, as our particle appears to 'scatter backwards in time ' when it encounters a 
lump of space-time curvature". However, when the heuristic argument given above 
holds, this situation is actually reversed: The mere co-existence of Q and ( | ) phys im- 
plies the existence of a unique preferred decomposition, and there is a well-defined 
notion of a conserved (not necessarily positive) probability-type density, provided 
by the scalar product Q (namely the integrand of Q2.8p ). What is missing in this 
framework is of course the notion of causality underlying the decomposition. Maybe 
one could say that the group average applied in this quantization procedure has 
singled out a decomposition by some kind of average of global type. 

In the case of the flat Klein-Gordon equation we actually know the above argu- 
ment to hold, and in addition we know that K 2 is a multiple of the identity operator 
(( | ) phys therefore being normalized such that K 2 = 1). This implies that ( | ) phys , 
when restricted to IH ± , is just ±Q. We do not know whether this feature carries 
over to more general cases. In general, we can only give a symbolic expression for K. 
As already mentioned in Section 2, physical states are represented as ^ phys = 8(C)ij) 
with C the wave operator in ( |2.2|) and ip a more or less arbitrary function (not satis- 
fying the wave equation). Defining an auxiliary operator K by {Kip) hys = 5(C)Kip, 
one finds 

K = - l -{6^n a V Q }6(C), (8.2) 

where E is a hypersurface on which Q is computed, S^y) the (^-function on E, 
such that f M dfj,SY:^p = f s dT*ip for any (test) function <p on A4, and dT> the scalar 
hypersurface element on S. The unit normal to S is n a , so that c?E Q = <i£n Q is 
the hypersurface element as used in ( |2.8|) . The expression ( |8.2| ) — when treated 
rigorously — should be independent of S. In flat space, its Fourier transform may 
be written down explicitly, exhibiting the expression 5{k^k^ + m 2 ). This is identical 
to | k$ l (5(kQ — oj£) — 5(k + ^)) with ujj: = (Ar + m 2 ) 1//2 , which illustrates the origin 
of the change of sign in the negative frequency sector. 
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As stated above, it is unclear how our heuristic argument leading to K relates 
to the formalism poposed in the main part of this paper, and how an analyticity 
condition might appear. The main condition for the refined algebraic quantiza- 
tion scheme to work is the self-adjointness of the wave operator C with respect to 
the scalar product (|2.7| ), and this is essentially a condition on the metric and the 
potential, without any obvious relation to analyticity issues. 



9 Outlook 

There is another important issue that seems worth being pursued. In this paper we 
have confined ourselves to the case that the potential U is positive everywhere in 
minisuperspace. As a consequence, the classical trajectories (whose tangent vector 
squared has the opposite sign as U) were always timelike, and we could clearly 
distinguish between incoming and outgoing modes at the classical level. Technically, 
we have used \/U in various expressions. However, in most realistic models the 
potential may take both signs. It is thus desirable to know how far the findings of 
this paper may be generalized for these models. One might still have a classically 
well-defined notion of outgoing and incoming traejctories (as long as M. is foliated by 
spacelike hypersurfaces) , but the classical trajectories would be spacelike in the U < 
domains (which are sometimes called "Euclidean"). Thus, a classical trajectory 
could start with outgoing orientation, dive into a Euclidean domain and re-enter the 
U > domain with incoming orientation. This might be a hint for the breakdown 
of the framework developed in this paper. In particular, evolution equations like 
( |3.23| ) would not easily be defined, and one might have to invoke the mathematics of 
signature changing metrics. Also, one could try to save the notion of decomposition 
by making it a local concept. On the other hand, the hints at analyticity playing a 
deep role in our framework open another possibility. Since analytic functions exist 
"as a whole", one could think of applying our formalism to the U > region and 
look whether it naturally "extends" to the whole of Ai (maybe in some complexified 
sense). Note that, for example, by describing a hypersurface E$ in terms of some 
(real) analytic function in the domain U > 0, a (complex) global behaviour is fixed 
at the same time. This might help clarifying the role of Euclidean domains and the 
complexified concepts (like Euclidean trajectories and Euclidean action) that are 
sometimes used in this field. 

As a last issue we mention that the most ambitious goal would of course be to 
carry over our proposal to the framework of the Wheeler-DeWitt equation in full 
(rather than mini-) superspace. In the metric representation p4| , one would have 



to deal with functional differential equations, and the analogues of evolution equa- 
tions like ( |3.23|) and (|3.38|) would be of the Schwinger-Tomonaga rather than the 



Schrodinger type, i.e. involve the feature of many-fingered time. Also, regulariza- 
tion issues are likely do make this a very hard job. As an alternative, one could try 
to work in the connection representation |25| and ask whether the notions of Klein- 
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Gordon type scalar product and decomposition can be given any sense therein. 



Acknowledgments 

I would like to thank Donald Marolf and Jose Mourao for stimulating my interest 
in the refined algebraic quantization program during the early stages of this work. 
Helmut Rumpf pointed out to me that in certain cases the operator K defined in 
Section 8 has in its point spectrum |T4]j . Also, I am indebted to Helmut Rumpf and 
Helmuth Urbantke for bringing the predecessors of the refined algebraic quantization 
program fTTJ-JTl] to my attention. 



References 

[1] J. J. Halliwell, "Introductory lectures on quantum cosmology", in: S. Cole- 
man et. al. (eds.), Quantum cosmology and baby universes, World Scientific 
(Singapore, 1991), p. 159. 

[2] S. W. Hawking and D. N. Page, "Operator ordering and the flatness of the 
universe", Nucl. Phys. B 264, 185 (1986). 

[3] A. Ashtekar, J. Lewandowski, D. Marolf, J. Mourao and T. Thiemann, "Quan- 
tization of diffeomorphism invariant theories of connections with local degrees 
of freedom", J. Math. Phys. 36, 6456 (1995). 

[4] D. Bohm, "A suggested Interpretation of the Quantum Theory in Terms of 
"Hidden" Variables. I" and "II", Phys. Rev 85, 166 and 180, (1952); J. C. Vink, 
" Quantum potential interpretation of the wave function of the universe" , Nucl. 
Phys. B 369, 707 (1992); Y. Shtanov, "On pilot wave quantum comsmology" 
(earlier title: "Pilot wave quantum gravity"), preprint |gr-qc / 9503005| ; A. Blaut 



and J. Kowalski-Glikman, "Quantum potential approach to class of quantum 
cosmological models", Class. Quantum Grav. 13, 39 (1996). 

[5] A. Vilenkin, "Interpretation of the wave function of the universe", Phys. Rev. 
1) 39,1116 (1989). 

[6] A. Higuchi, "Linearized quantum gravity in flat space with toroidal topology", 
Class. Quantum Grav. 8, 2023 (1991). 

[7] D. Marolf, "Quantum Observables and Recollapsing Dynamics", Class. Quan- 
tum Grav. 12, 1199 (1995). 

[8] K. Kuchar, "Canonical methods of quantization", in: C. J. Isham, R. Penrose 
and D. W. Sciama (eds.), Quantum Gravity 2. A Second Oxford Symposium, 
Clarendon Press (Oxford, 1981). 



42 



[9] A. Higuchi, "Quantum linearization instabilities of de Sitter spacetime: II", 
Class. Quantum Grav. 8, 1983 (1991); "Quantum Linearization Instabilities of 
de Sitter Spacetime", in: B. L. Hu and T. A. Jacobson (eds.), Directions in 
General Relativity, Vol. 2, Cambridge University Press (Cambridge, 1993), p. 
146. 

D. Marolf, "Refined algebraic quantization: Systems with a single constraint", 
preprint |gr-qc/9508015| , to appear in Banach Center Publications. 

O. Nachtmann, "Dynamische Stabilitat im de-Sitter-Raum" , Sitz. Ber. Ost.Ak. 
d. Wiss. 7/176, 363 (1968); "Continuous Creation in a Closed World Model", 
Z. Physik 208, 113 (1968). 

H. Rumpf and H. Urbantke, "Covariant 'In-Out' Formalism for Creation by 
External Fields", Ann. Phys. (N. Y.) 114, 332 (1978). 

H. Rumpf, "Mass-analytic quantization, uniform acceleration, and black-hole 
space-times", Phys. Rev. D 28, 2946 (1983) and references contained therein. 

H. Rumpf, "Self-adjointness-based quantum field theory in de Sitter and anti-de 
Sitter space-time", Phys. Rev. D 24, 275 (1981). 

C. Kiefer, "The semiclassical approximation to quantum gravity", in: J. Ehlers, 
H. Friedrich (eds.), Canonical Gravity: From Classical to Quantum, Springer 
(Berlin, 1994), p. 170. 

C. Bertoni, F. Finelli and G. Venturi, "The Born-Oppenheimer Approach to 
the Matter-Gravity System and Unitarity", preprint |gr-qc / 96040 1 1| . 

K. Kuchaf, "General relativity: Dynamics without symmetry", J. Math. Phys. 
22, 2640 (1981). 

R. M. Wald, "Proposal for solving the 'problem of time' in canonical quantum 
gravity", Phys. Rev. D 48, R2377 (1993); A. Higuchi and R. M. Wald, "Ap- 
plications of a new proposal for solving the 'problem of time' to some simple 
quantum cosmological models", Phys. Rev. D 51, 544 (1995). 

F. Embacher, "Quantum cosmology in the energy representation", preprint 
UWThPh- 1996-32, also |gr-qc/960502Tl , to appear in Nucl. Phys. B. 

F. Embacher, "On the interpretation of quantum cosmology", preprint 
UWThPh- 1996-22, also |gr-qc/96050i9| . 

N. D. Birrell and P. C. W. Davies, Quantum fields in curved space, Cambridge 
University Press (Cambridge, 1982). 

M. Abramowitz and I. A. Stegun (eds.), Handbook of Mathematical Functions, 
Dover Publications (New York, 1972). 



43 



[23] D. Marolf, private communication. 

[24] J. A. Wheeler, "Superspace and the nature of quantum geometrodynamics" , 
in: C. DeWitt and J. A. Wheeler (eds.), Batelles Rencontres, Benjamin (New 
York, 1968), p. 242; B. DeWitt, "Quantum Theory of Gravity. I. The Canonical 
Theory", Phys. Rev. 160, 1113 (1967). 

[25] A. Ashtekar, Lectures on Non-perturbative Canonical Gravity, World Scientific 
(Singapore, 1991). 



44 



